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METHOD FOR DETERMINING INDIVIDUAL HYDROCAR- 
BONS IN MIXTURES OF HYDROCARBONS BY MEAS- 
UREMENT OF FREEZING POINTS ?? 
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ABSTRACT 


z A method is described for determining individual hydrocarbons in mixtures of 
“MM hydrocarbons by the measurement of freezing points. Experimental data are 
given for the determination of the four Cs aromatic hydrocarbons. The uncer- 
tainty of the determination of each hydrocarbon is near +1 percent of the total 
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I. INTRODUCTION 


Street _— —_ 


This paper describes a general method for determining individual 
hydrocarbons in mixtures of hydrocarbons by measurement of freezing 
points. The method is based on the fact that the amount of a 
hydrocarbon H in an unknown mixture of hydrocarbons can be de- 
termined by measurement of the lowering of the freezing point of a 
sample of pure H caused by the addition to it of a given amount of 
the unknown mixture. If the unknown mixture contains none of the 
hydrocarbon H, the lowering of the freezing point will have its greatest 
value, which will be substantially that called for by the ideal or dilute 
solution laws; whereas, if the unknown is composed entirely of the 
hydrocarbon H, the lowering of the freezing point will be zero. 

In the exact application of the method, it is important to show that 
the solution of the hydrocarbon H with the unknown is, within the 
desired limits, ideal or sufficiently dilute, so that the ideal law of the 
lowering of the freezing point may be applied. If the solution of the 
hydrocarbon H with the unknown is not ideal or sufficiently dilute, 
it will be necessary to separate the unknown mixture beforehand by 
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an appropriate fractionating process into two or more portions, each 
of which will form substantially ideal, or sufficiently dilute, solutions 
with the appropriate pure hydrocarbons. 


II PREPARATION OF THE SAMPLE 


In order to utilize the freezing point observations in the desired 
manner, it is important that each mixture of hydrocarbons examined 
be one that forms an ideal or sufficiently dilute solution with each of 
the hydrocarbons to be determined. In general, this requires that the 
unknown mixture of hydrocarbons be composed preponderantly of 
hydrocarbons of a similar or nearly similar type, as for example, (a) ' 
a mixture of paraffins plus cycloparaffins or of either of these separately. 
or (b) a mixture of olefins, or (c) a mixture of aromatics. Such mix- 
tures will in general be ones most likely to follow, within the desired 
limits, the ideal law of the lowering of the freezing point. On the other 
hand, mixtures containing hydrocarbons of different type (considering 
for this purpose paraffins and cycloparaffins to be of similar type) may 
be expected to show significant departures from the behavior of an 
ideal solution. Such mixtures would require appropriate fractionation 
by type before proceeding with measurements of the lowering of the 
freezing point. 

The method of analysis described in this paper places no restriction 
or limitation on the boiling range of the material to be examined, 
although obviously the precision is likely to be better for a narrow- 
boiling fraction than for a very wide-boiling fraction, such as the entire 
gasoline fraction. 

The amount of hydrocarbon material required is, for each determina- 
tion, about 40 ml of the “pure” hydrocarbon and 2 to 3 ml of the 
unknown. 

III. METHOD 


Given an unknown mixture of hydrocarbons of similar type, it is 
required to determine the amount of a given hydrocarbon H which 
is in the unknown mixture. Measurement is made of the freezing 
point of a sample of ‘‘pure”’ H, and also the freezing point of a dilute 
solution of the unknown mixture in ‘pure’ H. The molecular 
weight of the unknown is determined, if necessary, by measuring the 
lowering of the freezing point of another hydrocarbon, of the same 
type as the unknown mixture but which is not a component of the 
unknown mixture, on the addition to this other hydrocarbon of a 
measured small amount of the unknown. Let 


N?=the mole fraction of the hydrocarbon H in the given lot of 
“pure” H. (This will be slightly less than unity, in general.) 

N*=the mole fraction of the hydrocarbon H in the solution made 
by adding a small amount of the unknown to “pure” H. 

N*“=the mole fraction of the hydrocarbon H in the unknown 
mixture. 

m?=the mass of the sample of ‘pure’? H used in making the 
solution of ‘‘pure” H and the unknown. 

m*“=the mass of the unknown used in making the solution of 
“pure” H and the unknown. 

M’=the molecular weight of ‘‘pure’’ H, taken to be that given 
by the molecular formula of hydrocarbon H. 
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M*“=the molecular weight of the unknown, determined, if neces- 
ren from a lowering of the freezing point as mentioned 
above. 

t,=the freezing point of hydrocarbon H for zero impurity, in 
degrees centigrade. 

t?=the freezing point of the given lot of “pure” H, which, in 
general, will be slightly pews than t,. 

t=the freezing point of the solution of “pure” H with the 
unknown. 


It follows that 


Pp Ue Dp u 
pt NM ( anti) (1) 
which is equivalent to 
?\ (/ M* 
Nan—(7 ae )N?—N?)- (2) 
The values of N? and NV’ are given by the following relations [1, 2, 3): 
—In N?=A(t,—t/)[1+Bt,—t/)] (3) 
—hn M=A(t,—t/)[1+Bty—t,)I, (4) 
where 
A=AH,/RT;,,’, (5) 
B=1/T,,—AC,,/24Hq, (6) 
T, =the freezing point in degrees Kelvin of pure H with zero 
impurity, 


R=the gas constant per mole, 
4H,,=the heat of fusion of pure H at T',, 


AC,,= the heat capacity of pure H in the liquid state less the heat 
capacity of pure H in the solid state, at 7’. 


In eq 2, which gives explicitly the value of the desired quantity N*, 
the values of N? and NV’ are computed from eq 3 and 4, respectively, 
using the measured values of t? and ¢t,, with the given value of 
ty. It should be noted that for high precision the difference between 
t? and ¢t,? must be as precise as possible, whereas the value selected 
for ty, is of secondary importance, since it is used primarily as a 
reference point. In fact, if the sample H is sufficiently pure with 
N? very nearly equal to unity, then eq 2, 3, and 4 may be simplified 
without significant error by substituting 1 for N? and t/? for ty, as 
follows: 


Non —-(Z)( 555) —N’) (2a) 
—In N°=0; N’=1 (Ba) 
—In N= A(t?—t)[1+B(t?—t/)]. (40) 
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In eq 2 (and 2a), the masses m? and m* are obtained from weighings 
on a balance, the molecular weight M? is taken as known and equ 
to that of the hydrocarbon H, and the molecular weight of the up. 
known mixture, M*, is determined as follows: 

Given a hydrocarbon G which is not a component of the unknow 
mixture but which forms with it a substantially ideal or sufficiently 
dilute solution. Measurement is made of the freezing point of a lo 
of “pure” G, which should contain not more than 0.005 mole fraction 
of impurity, and of the freezing point of a dilute solution of the 
unknown in G. Let 


t?” =the freezing point of the lot of “pure” G, 

t;* =the freezing point of a dilute solution of the unknown in G, 

m? =the mass of “pure” G in the solution, 

m“ =the mass of the unknown in the solution, 

M?’=ithe molecular weight of G, taken as known from its molec. 
ular formula, 

M*=the molecular weight of the unknown mixture. 


If the freezing point of “pure” G, t/, is not too much below the freez- 
ing point of G for zero impurity, t,, then the molecular weight of the 
unknown is obtained with sufficient accuracy from the following 
equations: 


—In N= A(t?—t/)[1+B(i?—t,)] 
and 
M*=M?(m*/m?)N*/(1—N*). 


In eq 7, the constants A and B are the same as before, except that 
they now apply to pure hydrocarbon G. 

he possible precision of the method in the case of the four (, 
aromatics is as follows. It is necessary to determine the difference 
between t? and t,; within 0.01° C, if the amount of the given com- 
pound is to be determined to within 0.005 mole fraction, when the 
relative amounts of the “pure” hydrocarbon and the unknown are 
selected so as to give a lowering of about 2° C when none of the given 
compound is present in the unknown. 


IV. EXPERIMENTAL PROCEDURE AND RESULTS 


The use of this method is here illustrated by application to the four 
C, aromatic hydrocarbons. Since it has been shown [8, 9]* that the 
four C, aromatic hydrocarbons follow the ideal solution laws as regards 
the liquid-solid equilibrium, it appeared desirable, for this system, to 
use as known test mixtures simply the four pure components in turn, 
determining in each case a component different from the one being 
examined so that the known would contain none of the component 
being determined. In these experiments, the lowering of the freezing 
point would, for the given molal ratios used, be a maximum, and 
departures from the theoretical values would be more easily recognized. 

In the experimental part of the work, the determination of freezing 
points was made with the apparatus and procedure described in 
reference [3], with the following additions to the procedure: 


4 Figures in brackets indicate the literature references at the end of this paper. 
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(a) For best results in the case of solutions, it was found necessary 
to extend the observations on the equilibrium portion of the time- 
temperature freezing or melting curve to cover the largest possible 
fraction of the material crystallized or melted. 

(b) Instead of a visual extension of the equilibrum portion of the 
curve to “‘zero” time on freezing (the time at which crystallization 
would have begun in the absence of undercooling), or to the corres- 
ponding time on melting, the geometrical method described in a 
recent report of the' API Research Project 44 at this Bureau was 
used [4]. : 

For certain hydrocarbons, it is found that the time-temperature 
freezing curve obtained in the manner referred to has no portion 
recognizable as representing thermodynamic equilibrium between 
solid _ liquid. In such cases, the time-temperature melting curve 
is used. 

To test the over-all accuracy of the method, determinations were 
made of the amount of ethylbenzene, of p-xylene, and of m-xylene 
in “pure” o-xylene (0.9940 mole fraction pure), and the amount ef 
o-xylene in ‘‘pure” p-xylene (0.9948 mole fraction pure). The 
amounts determined (see table 1) were, respectively, 0.0090,—0.0012, 
0.0138, and 0.0034 mole fraction, which values are to be compared 
with zero as the theoretical value. 

A convenient way of tabulating the experimental results and 
calculations is illustrated in table 1, where the columns give, in order, 
the following quantities wherever appropriate (see footnotes to the 
table for further explanation): the freezing point for zero impurity, 
N,=0; the resistance in ohms of the given platinum resistance ther- 
mometer at the given freezing point; the difference in resistance of the 
platinum resistance thermometer at the given freezing point and at 
the freezing point for zero impurity; the change in resistance with 
temperature for the given therniometer at the given mean tempera- 
ture; the value of ty,—t,, obtained as the quotient (R,;,—R,) /(dR/dt) ; 
the value of A as defined by eq 5; the value of B as defined by eq 6; 
the ratio of the masses of the ‘‘pure’’ substance and the unknown; 
the ratio of the molecular weights of the unknown and the “pure” 
substance; and the mole fraction of the given hydrocarbon in the 
“pure” substance, in the solution, and in the unknown. Table 1 also 
gives the values of A and B, as defined by eq 5 and 6, respectively, 
for benzene. 

For these small differences, the calculation of the value of the 
lowering of the freezing point may be made most conveniently and 
precisely by taking the difference in the two values of the resistance, 
R, of the resistance thermometer, and dividing by the value of 
dR/dt for the mean temperature. (If a thermoelement is used, the 
difference in the two values of the emf, E, is divided by the appro- 
priate value of dE#/dt.) For a platinum resistance thermometer for 
which the constants Ro, 6, 8, and c (=(Rio—Rp)/Ro) are given, 


dR/dt=Rge[(1+0.015) — (2 10~*)5t— (48 10-8) (t—75) #]. 


The foregoing equation is applicable for temperatures below 0° C. 
For temperatures above 0° C, the equation is simplified to 


dR/dt=R,e[(1+0.018) — (2 10) st]. 


582610—44-2 
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With these equations, values of dR/dt for unit values of t covering 
the range of interest may be tabulated and linear interpolation used 
for obtaining values for intermediate temperatures. 

Application of the method has also been made [9] to the determina. 
tion of the amounts of each of the four Cs, aromatic hydrocarbons jy 
five different unknown mixtures consisting of C;, Cs, and Cy aromatic 
hydrocarbons, and the amounts of ethylbenzene and p-xylene in 
seven other different unknown mixtures consisting of C;, Cy, and(, 
aromatics. For each of the first five of these unknowns, the following 
independent data were available to serve for comparison with the 
results obtained by the method of measuring freezing points: 

(a) From analytical distillation data the relative amounts of the 
C,, Cs, and C, aromatics (from which values the mean molecular 
weights could be calculated, as well as the percentage of the sum of 
the four C, aromatics). 

(b) From analytical distillation data, the percentage of o-xylene in 
each of the unknowns. 

The results obtained on one of the unknowns by the present method 
are given in table 1. 


TABLE 1.—Ezperimental data on the determination of the four Cs 
aromatics in an unknown miature of aromatics 




































































Material ty, | Ry | Ry—Ry | dR/dt tyy—ty 
ETHYLBENZENE 
°C Ohms Ohms _ | Ohms/degree °¢ 
No=0 ®_____- ny tee , ’ — 94. 950 15. 7212 | eee D 0.000 
“Per Fi oi Scabbbbbasctoe RE epee am & Pet | ----| 15.7176 . 0036 0. 10509 - 034 
“Pure’’+-o0-xylene__....---.--- Seid aiitacainhin 15. 5094 . 2118 . 10514 2.014 
“Pure”’+unknown.-_-_-__-.-- oie Bee . 1853 . 10514 1, 782 
p-XYLENE 
No=0 ®__..-.__._- a c San 13. 260 26. 8621 OB bas A 0.00 
td, eS eee rae teh ceside assy pete lebp<esaeay 7m eee . 0210 0. 10116 P.) 
“‘Pure’’+-o-xylene__. aap fala # PED, secstln ta « 26. 5995 . 2626 . 10120 2. 505 
“Pare aneeowe 20 ki Sede als ds Le ee . 2074 - 10119 2.050 
M-XYLENE 
og a EE Pe OP Me SS SS Bay oe 55 ES Reena —47. 890 20. 6183 | OOOO {iis .ceccace 0.00 
a EE Ee Ree ne enn Bie, Sree 20. 4836 . 1347 0, 10314 1, 306 
Sc ty a, RES R IRS SERS APOE ESSA Nek ire 20. 2735 . 3448 . 10318 3. 342 
“Pape? SOIR 5 ok 5. es edab bcs. nck sc. scodd ekphnceac 20. 3611 . 2572 . 10316 2. 498 
| | 
0-XYLENE 
ROD Os ihe red ebhec dp thd - cancetrinwspwane shops — 25. 190 22. 9503 0. Py eer 0, 000 
SNE cite ee ethos san nh +8) rere atietn ial ngies isin im . 0232 0. 10234 ami 
“Pure”’+p-xylene.........-..-- TALE vith se Jods dadbtd 22. 7149 . 2354 . 10238 2. 20 
“Pure’+unknown. .......--- rn reste SPREE 22. 7489 . 2014 . 10237 1, 961 
| a 
BENZENE 
C2 ee RNS Wl Se A ER EA ELAS ee m= ERA US SPE. 
Be cen aed | ulna A | | 








See footnotes at end of table. 
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TaBLE 1.—Ezperimental data on the determination of the four Cs 
aromatics in an unknown mizture of aromatics—Continued 










































































Material A B m>/ms | M+*/M» Np | Ne | N« 
ear | | 
ETHYLBENZENE 
| Mole Mole Mole 
Degree! | Degree-' | fraction | fraction | fraction 

ih, sonicedinineteriieitts 0.03454 | 0.00424 |......-..-]......---. TS i TRE ES SeenON Ie 
OMT. g no dnchp ner ddbuaxkhirabdosedsnlbgueusen’slantcuxEaeelercesncc ad |. RA fees . 
“pare GMMR 5. wo ol Sop ek. 243 Pek kk 13. 871 1. 0000 .9988 | 0.9322 0. 0090 
Beas AES” REI & AS RRR 14.072 | 1.0085 . 9988 . 9405 . 1136 

p-XYLENE 
Pie ss Somes ae ' 
Nees... acc ccesctesacesssis] OOM O.90878 | ce Cet det ieele. Bieas Th SRR 7 ania 
| ee RE EIS GN REC EN. weeeeerene|eeesenccan[oonzennne-| 0.0068 |.-.-------|.------- 3 
“Pure’+-o-xylene....---.------ ORES (ani Meee May yee .9948 | 0.9365 | —0.0012 
“pyre”+unknown.............|........-- SX # --| 17.236 | 1. 0085 | 9048 | . 9496 1641 

™M-XYLENE 
vv ae ARR ep epee tegs ve Rema OO) Ge ha TA ae Soe Be Re 
Fe th ae Sembee k Spann: syed rating Sia: yon a apOO ee CORDS i  *S ee Shee 
“Pyre”-+-o-zylene......--------|..-----.-- ae as 16.762 | 1.0000 .9644 | 0.9109 0. 0138 
“Pure”+unknown.......-....- ee | se BR | 16334} 1.0085 | [9644 | 9329 4142 
! i ! ! 

0-XYLENE 
See ee, ORM Ta aan 0.02660 | 0.00312 |..........]|_- ik va | 1. 0000 | ae # | eee 8 
5 TES ES AEE TOE CE BREN = Ee TEE eel CAN Ge RR EU Se enna 
“Pyre’+-p-xylene........--.---|------ 17. 442 1. 9940 | 0.9403} 0.0034 
“Pure” +unknown... facsoqewese Jenwene te- | 17.091 | 1.0085 | 9940 . 9487 | 1688 

BENZENE 
Weir ho eee es | 0. 01523 | 2 ER Se, Eanes | eo. eres Se | Eeaee! > | 














_* This row, labeled N2=0, gives values for the substance with zero impurity, that is, having Na=0, where N3 
is the sum of the mole fractions of all the impurities. The indicated freezing point for zero impurity is the 
— value, taken from the tables of the API Research Project 44 at the National Bureau of Standards. 
_ > This row, labeled ‘‘Pure’’, gives values for the material actually used for the ‘‘Pure”’ sample. Its purity 
isgiven in the column headed Np. 


The molecular weight of the unknown, determined from the lower- 
ing of the freezing point in benzene, is 107.31, which is to be compared 
with the value 107.06 determined from the data of the analytical 
distillation. In a check experiment, the molecular weight of a sam- 
ple of pure ethylbenzene was determined from the lowering of the 
freezing point in benzene to be 106.45, as compared with the theoreti- 
cal value of 106.16. The sum of the percentages of the four C, 
aromatics was determined to be 86.1 +2.4 percent by volume by 
measurement of freezing points, which value is to be compared with 
85.0 +0.5 percent by volume obtained from the data of the analytical 
distillation. The amount of o-xylene determined by measurement of 
freezing points is 16.9 +1.2 percent by volume, which is to be com- 
pared with the value 17.6 +1.0 percent by volume obtained from 
the data of the analytical distillation. 

Figure 1 shows a time-temperature freezing curve used to deter- 
mine the freezing point of a solution of the unknown in o-xylene. 
Figure 2 shows a time-temperature melting curve used to determine 
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the freezing point of a solution of the unknown in ethylbenzene 
Figure 3 shows a time-temperature freezing curve used to determing 
the freezing point of a solution of the unknown in benzene, for th 
determination of molecular weight. 


V. DISCUSSION 


From the results of the data so far obtained, it appears that, with 
the method and procedure described here, individual eromuelin hydro. 
carbons can be determined in a mixture of aromatic hydrocarbons 
with an uncertainty for each compound near +1 percent of the 
sample. 
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Individual Hydrocarbons in Mixtures of Hydrocarbon 


Previous measurements [3] have shown that mixtures of paraffins 

and cycloparaffins in dilute solution follow the ideal solution law of 
the lowering of the freezing point. It follows, therefore, that the 
present method can be used to determine individual paraffin or cyclo- 
paraffin hydrocarbons in mixtures of paraffin and cycloparaffin hydro- 
arbons. 
Although no information is available with regard to olefin hydro- 
carbons and other types, there appears to be no reason why the 
method can not be used to determine individual olefins in a mixture 
of olefins, etc. ; 

When the present investigation was nearing conclusion, references 
were found to the following three papers in the literature that men- 
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tion or describe a method of analysis similar in principle to tha 
described in the present paper: Norris and Rubinstein [5], Norris and 
Vaala [6], and Ibing [7]. The first two of these papers [5, 6] ar 
concerned with other problems and barely mention the method. 
In the third paper, [bing [7] describes his procedure and apparatus 
in detail, uses in his calculations only the simple limiting form of the 
law of the lowering of the freezing point, and does not define the 
freezing point nor tell how he determined it from his time-temperature 
observations. 
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TIME IN MINUTES 


Figure 3.—Time-temperature freezing curve for benzene plus an “unknown” 
mixture of aromatic hydrocarbons. 


The scale of ordinates gives the resistance of the platinum resistance thermometer in ohms. The scale of 
abscissas gives the time in minutes. 
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TURE FREEZING AND MELTING CURVES AS APPLIED 
TO HYDROCARBONS ?? 


By William J. Taylor * and Frederick D. Rossini 


ABSTRACT 


A method is described for determining analytically or graphically, from appro- 
priate time-temperature freezing and melting curves obtained on hydrocarbons, 
the freezing point of a given substance, and within certain wide limits, the freezing 
point of that substance for zero impurity (liquid-soluble, solid-insoluble). 
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I. INTRODUCTION 


In a recent paper [1] * dealing with time-temperature freezing and 
melting curves obtained on hydrocarbons, it was shown how, when 
significant undercooling occurs on freezing, the observations can be 
analyzed. to yield, in proper relation to the subsequent equilibrium 
part of the curve, the “zero” time, that is, the time at which crystalliza- 
tion would have begun in’the absence of undercooling. In that paper, 
visual extrapolation of the equilibrium part of the curve, to ‘zero’ 
time on freezing and to the corresponding time on melting, was used 
to determine the freezing point, defined as the temperature at which 


' This investigation is part of the work on Research Project 44 of the American Petroleum Institute, from 
whose research fund financial assistance has been received. 

‘Presented before the Division of Petroleum Chemistry at the meeting of the American Chemical Society 
st Cleveland, Ohio, April 5, 1944. : 

‘Research Associate at the National Bureau of Standards, representing the American Petroleum Institute. 

‘ Figures in brackets indicate the literature references at the end of this paper. 
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an infinitesimal amount of crystals of the major component is jp 
equilibrium with the liquid.° 

The present paper describes a method, which may be used ana- 
lytically or graphically, for extrapolating the equilibrium part of , 
time-temperature freezing or melting curve to give the freezing point, 
T,, of the given hydrocarbon substance. There is also included 
description of a method for deducing, within certain wide limits, from 
the same observations on the given substance if carried over a suff- 
ciently large fraction of material crystallized or melted, the freezing 
point for zero impurity, T,,. When the heat of fusion is known or can 


be estimated, the value of 7;,—T; serves to give the amount of 
liquid-soluble, solid-insoluble impurity in the given substance. 


II. DERIVATION OF THE BASIC EQUATIONS 


The definitions, abbreviations, and symbols used in this paper will 
be exactly the same as those of the previous paper [1], except that the 
subscript f attached to the time z, as in z,, denotes the time corre- 
sponding to the freezing point, 7, Some of the previous definitions 
and equations are repeated below for the sake of clarity and emphasis. 

For an ideal, or a sufficiently dilute solution, the thermodynamic 
relation between the temperature of equilibrium and the composition 
of the liquid phase, for the equilibrium between a liquid phase consist- 
ing of a major component and one or more other components and a 
crystalline phase of the major component alone, is [2, 3] 


—In Nj=—In (1—N,)=Ni(1+1/2Nit+ . » = 
(AH,,/RT';,7) (T,,.—T) {1+ (1/T,—AC,,/24H,) (7,,—T) + oe be (1) 


In eq 1, the symbols are defined as follows: 


N,=the mole fraction of the major component in the liquid phase; 
N.=(1—N,)=the sum of the mole fractions of all other com- 
ponents in the liquid phase; 
R=the gas constant, per mole; 
T;.=the absolute temperature of the freezing point of the major 
component when pure (that is, when N,=1 or N,=0); 

A4H,,=the heat of fusion, per mole, of the major component in the 
pure state at the temperature 7',,; 

AC,,= the heat capacity per mole, of the pure liquid less that of the 
pure solid, for the major component in the pure state at 
the temperature 7',,; and 

T=the given temperature of equilibrium. 


The assumptions involved in the derivation of eq 1 are the following: 

(a) The liquid solution is ideal, or sufficiently dilute, so that the 
fugacity of the major component in the liquid phase is proportional 
to its mole fraction. 

(b) The heat of dilution is negligible over the range of concentration 
from N,=1 to the given value of N,. 


5 It is important to note that this temperature, when properly measured, must have the same value whether 
approached from a higher temperature as in cooling or from a lower temperature as in warming. 


Freezing and Melting Curves of Hydrocarbons 199 



















(c) The heat of fusion of the major component in the pure state at 
the temperature 7’ may be represented by a power series expansion, as 


4 AH,=AH,,+AC,,(T—Tp)1+ .. |. (2) 
4 Introducing the following additional abbreviations, 

1 A=AH,,/RT,,?, and B=1/T;,— AC,,/24H 4, (3) 
* expanding —In (1—N), and rearranging, eq 1 becomes 


Ty—T=U+12N+ .. -BTy— T+ ob @ 


Neglecting terms in N° and higher, eq 4 becomes 


Ty—-T=SU+1/2Ne+ . \-[B *) es: A (5) 
or 


In a system of the type under investigation, in which the solid phase 
is restricted to the major component, let n,* represent the total num- 
ber of moles of the major component (in both the liquid and solid 
phases), me represent the total number of moles of the other compb- 
nents (in the liquid phase), and n, represent the number of moles of 
the major component in the liquid phase under given conditions. 
Then the number of moles of the major component in the crystalline 
phase under the same conditions is n,*—n,. Let r denote the fraction 
crystallized of the total number of moles of all components, that is, 


r= (n*—n)/(m*+n). (7) 


The mole fraction, in the liquid phase, of the sum of all components 
other than the major components is defined by 


Na=MNa/ (+72). (8) 


Let N.* represent the value of N, when no crystalline phase is present, 
that is, when m=7,*. Then 


N2* =N2/(m* +12). (9) 


— of eq 7, 8 , and 9 yields the following equation relating 
N, and r: 


N.=N,*/(1—7). (10) 
Letting 
a=N,*/A, and b=(1/2—B/A) (N2*)?/A, (11) 


and combining eq 10 and 6, there is obtained 


oh arial rn ae aie: =Ty—;*,[ 1+ 75+ a ‘b (12) 
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It is possible to further simplify eq 12 by using the general relatio, 


that, when 
b 2 
(Hist 


(abate) Slt - | (13) 


Assuming that the unwritten terms in eq 12 and 13 are negligible, or 
may be substituted one for the other without significant error, they 
eq 12 may be written 





In general, the fractional error introduced by making the approxi- 
mations leading to eq 14 will be of the order of 1/12 to 1/3 N,?, de. 
ending upon the value of B/A (see eq 3 and 4). That is to say 
ractional errors in N2, or correspondingly in 1;,,—T, will be about 


1/3 to 1 1/3 percent when N.=0.2, and about 1/12 to 1/3 percent 
when N,=0.1. Equation 14 therefore would appear to be practically 
useful only when the amount of solute or impurity is below several 
tenths in mole fraction. 

The next step is to express r, the fraction crystallized or melted 
in terms of the time. This is done by assuming that the experiment 
is one in which the rate of crystallization or melting of the major 
component is constant with time.*® 
Denoting the time by the symbol z, and letting z=z, when r=0, one 
has the relation 

r=k(z—2,), (15) 


where k is a constant characteristic of the experiment, positive in 
i for a freezing experiment and negative for a melting experiment. 
The time at which freezing begins, or melting is complete, is z,, pro- 
vided thermodynamic equilibrium is maintained throughout the 
experiment. The determination of z, for an actual experiment has 
been discussed in reference [1]. 

Combination of eq 14 and 15 yields the following relation between 
the time and the temperature of equilibrium: 


, 


a ’ 
Lets ~T—VG—2,)' (16) 


where a’=a/(1—6/a) and k’=k/(1—}/a). 


Equation 16 relating 7 and z and eq 14 relating 7’ and r are the 
equations of a rectangular hyperbola. The asymptote obtained by 
setting z=— © in eq 16 is T=T,,. F 
Since there are three constants in eq 16, it is clear that three points 
on a portion of the curve satisfying eq 16 are sufficient to determine 
the values of the three constants, 7’, a’, and k’. However, if it is 


¢ This is essentially the case in experiments performed with the apparatus described in reference [1] whe0 
the range of freezing (or melting) temperature is sufficiently small in comparison with the head 0! 
temperature. 
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desired to test the accuracy with which the data satisfies the equation, 
a greater number of points are necessary. 

The following theorem concerning a rectangular hyperbola will be 
required. Consider any four points (7;, 2), (T2, 22), (Ts, 2s) and 
T,, 2), on the hyperbola defined by eq 16. There are three ways 
in which two chords can connect these points, each point being the 
terminus of only one chord. These are illustrated in figure 1, in 
which my denotes the slope of the chord connecting points (71, 2) 
and (72, 22), and similarly for m3, m3, and my. It is to be shown that 
the products of the slopes of the two chords is the same in each of 


the three cases 











M2734 = M1 3Meo4 = M4 4M 93. (17) 
Using eq 16, 
ES ah). —a’k’ 
mu=\ 2, ) [1k a2) —F 2] (18) 
and 
4)./\2 
M2M34 = (a k ) 


—¥@—2)l—F een —F @— ald — Pee 


Equation 17 follows immediately, since the expression on the right- 
hand side of eq 18a is unchanged in value by any interchange of the 
subscripts 1, 2, 3, and 4. 


Temperoture —— 








Time —» 


Ficure 1.—Construction to show three ways of connecting four points on a curve 
with two chords, with each point being the terminus of only one chord. 


III. DETERMINATION OF THE FREEZING POINT 
1. ANALYTICAL 


The procedure for determining the freezing point of the given 
substance analytically from the equilibrium portion of the time- 
temperature curve is illustrated by means of figure 2, Points G, H, 
and J are on the equilibrium part of the time-temperature curve, as 
defined by eq 16. The zero time, at which crystallization would have 
begun in the absence of undercooling, has been determined [1]. It is 
tequired to determine on the ordinate through the z ro time, z,, the 
freezing point 7’. 
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Figure 2.—Construction to show the analytical determination of the freezing point. 


From eq 17, it follows that 


GE Ln— “) T,—T. \- E ="), (19 
2g Z2y 217 Zh Zi 2 Znh— 2g 
Rearranging and using the relation 

(T,—T) =(T;—T,) + (T,—T)), (20 


one obtains 


~ 


FER) -(BERNGERNGED) 


Denoting by wu, v, and w the dimensionless ratios 


7 
OY at? 4 
. 4 ig 4 a 
r.=¥, 
Chaz IQ) 
<i nh 
—- 
w= a 21", (24 
eg 24) 


there is obtained, on substitution into eq 21, the relation 
T,—T;\ 


" (25) 
uvw—1 


T= T,+ 
Equation 25 gives the value of 7;, when there are known the tempert- 
ture and time corresponding to three points (as G, H, and J) on the 
equilibrium curve (as defined by eq 14) and the “‘zero”’ time, the time 





ont, 


(19) 


(25) 


eT a- 
the 
ime 
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at which crystallization would have begun in the absence of under- 
cooling.’ 

It is nearly always possible to select the point H equidistant in time 
between G and J, so that »=1. Furthermore, it will usually be pos- 
sible to space the points @, H, and J so that w is a simple ratio. Thus, 
if all three time intervals are ‘equal 


(2,—2y) = (Zn—2,) = (2:—2n) 


on 
o=8 26) 


T;=T,+ (T.—T,)/(8u—1) 


When the first time interval is twice as great as the other two intervals, 
then 


(2,—2y) =2 (Zn—2,) =2(21—2n) 


=? ! @7) 


T;=T,+ (T,—T;)/(u—1) 


It is obvious that a large value of w corresponds to a short extrapola- 
tion to 7, and vice versa. 

Although the curve in figure 2 has been drawn to simulate a freezing 
curve, the above equations apply equally well to a melting curve, since, 
for the same substance, the melting curve is the mirror image of the 
freezing curve when the time rate of melting is made the same as 
the time rate of crystallization. 


2. GEOMETRICAL 


The manner of determining geometrically the freezing point, 7, 
is shown in figure 3. 

There are given the “zero” time, z,, and the time and temperature 
corresponding to the three points, G, H, and J, on the equilibrium 
part of the time-temperature curve defined by eq 16. It is required 
to locate the freezing point, 7,, on the ordinate at the time z,. 

The solution is as follows: 

Draw AB through J parallel to the time axis. 

Draw AC perpendicular to AB, at the time 2z,. 
_ Draw a line through points @ and H, intersecting AB at E and 
AC at D. 

Draw a line through points J and H, intersecting AC at J. 

Draw JK parallel to DE. 


Draw a line through points K and G, intersecting AC at F. 
F is the desired point, representing the freezing point, T;. 





"It is worthwhile remarking that eq 25 is not limited to the determination of the freezing point = 
zin eq, 24 may be any given time, and eq. 25 then yields the corresponding temperature, 7’. It has 
implied in the foregoing discussion that the cooling curve before crystallization begins (or the wartning 
curve after melting is complete) is sufficiently steep so that zy May be estimated exactly enough from a 
preliminary visual extrapolation to Ty. It is clear that if this is not the case, 7’; may still be obtained with 
any required accuracy by successive applications of eq 25, 
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The proof, illustrated in figure 4, is as follows: 
Draw FT. 
Draw NP through H, parallel to AB. 




















Then fe ws V 
NH/NP=ND/NJ (28 
But i 
Al/AK=NH/NP (2) 
Therefore, 
ND/NJ=Al1/AK. (30) 
C ! 
(T,) m4 
J 
D 
t 
i (t)) F 
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K as a a 
f Time —> 


Figure 3,—Geometrical construction for determining Ty and T jo. 
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But, by eq 17, 


(ND/NH)(AF/AI])=(AF/AR)(NJ/NB), (31) 
which gives 
ND/|NJ=A1/AK, (32) 


in agreement with eq 30. 


IV. DETERMINATION OF THE FREEZING POINT 
FOR ZERO IMPURITY 


The method described ere for determining the freezing point for 
zero impurity, T',, can be properly applied only to those experiments 
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FieurE 4.—Proof of geometrical construction for delermining Ty and Ty, 
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in which the time-temperature observations extend over a sufficient] 

large fraction (probably of the order of % to 4) of material crystallized, 
Furthermore, since that portion of the curve which is required for 
determining 7’,, is also the portion in which departures will most likely 
occur in the relations upon which eq 16 is based, it is seen that the 
determination of even an approximate value of T,, will frequently be 
impossible for many experiments. However, for completeness, the 
description of the procedure for evaluating T,, by this method is given 
here so that it may be available for application to those experiments 
which do conform to the requirements. 


1. ANALYTICAL 


An equation for the freezing point for zero impurity, T fo» may be 
arrived at as follows. Returning to eq 16, it will be noted that 7, 
as defined there, approaches 7; as z approaches —. This may be 
regarded simply as an analytical property of eq 16, so far as the present 
argument is concerned.® 

Since eq 25 applies to any point on the curve defined by eq 16, we 
may imagine z, and 7; in eq 24 and 25 to be replaced by z and 7, 
respectively, corresponding to any point on the curve. Then as 2 
approaches — ©, w approaches unity, and eq 25 becomes 


T= 


1) 


(33) 


Since w has been eliminated from eq 33, T f. does not depend on the 


location of z; As pointed out following eq 25, it is nearly always 
possible to select the point H equidistant in time between G and J, 
so that v=1. 

The precision of the determination of the freezing point for zero 
impurity rests heavily upon the value of the product (wv), which can 
be seen from eq 22 and 23 to be the ratio of the slope of the chord 
HI to the slope of the chord GH. 


2. GEOMETRICAL 


The manner of determining geometrically the freezing point for 
zero impurity, T',, is shown in figure 3. 


After having determined the freezing point, 7;, as indicated in the 
preceding section, the procedure is as follows:° 


Draw a line through G parallel to AB and a line through J parallel 
to ACU, the two lines intersecting at L. 


Draw a line through points K and L, intersecting AC at M. 
M is the desired point, representing the freezing point for zero- 
impurity, 7’. 


8 The extension of the time in eq 16 to — may be easily pictured by considering a corresponding melting 
curve, which may be prolonged indefinitely in time by the continued addition of the pure major component 
in the solid phase to the existing solid phase as melting proceeds. This procedure serves to produce 8 liquid 
phase in which the number of moles of the major component approaches infinity, the composition of the 
liquid phase approaches N;=1 or Ns=0, and the temperature approaches the freezing point for zero im- 
purity, Ty. 

* The method of obtaining T;, described here has the advantage that 7;, is obtained on the same vertical 
axis as Ty, so that the initial freezing point lowering (7/,— 7) is immediately apparent. However, It 


should be pointed out that the value of 7, obtained depends only on the points G, H, and J, and not at 
all on the time chosen for zy. 
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The proof is as follows: 
Draw GJ. 


Since point F' is on the curve, eq 22, 23, and 33, for Zr may be 
rewritten with the three points G, H, and I or by F, G, and J. 


1 =1/+(35 iy) (34) 


(35) 


, (a2) (36) 


Equation 34 may be put in the form 


ae) , 
=1'¢’ . (37) 
T,,—T; 


where a ni r 
“ ee 


but from figure 4, 


7, mn 
AM/FM=IL/RL. 
(7 T, / / 


Also, from eq 35 and 36 and figure 4, 
u'v’ =M¢1/Mre, (39) 


where mg, and mpg are the slopes of the chords GJ and FG, respectively. 
Combining eq 37, 38, and 39, 


IL|/mex=RL | myo. (40) 
But both sides of eq 40 are equal to GL, from figure 4. 


V. UNCERTAINTY OF THE DETERMINATION OF THE 
FREEZING POINT AND OF THE FREEZING POINT FOR 
ZERO IMPURITY 


An expression may be derived for the error of 7 or T';,, as calculated 
from eq 25 and 33, for any given experiment. For simplicity, one 
may assume that there is no error in 2,, 2, Za, and 2, but that a 
uniform random error, ¢, is to be attached to 7;, 7, and T; (fig. 2). 
| The error in Jy, oy, may then be calculated from the usual propaga- 
tion of error equation 


F-CY CH w 


Equation 41 assumes that the errors in T,, 7, and 7; are, on the 
average, uncorrelated, and gives the error for 7; in terms of the 
average of these errors, o. 

The derivatives in eq 41 may be obtained by differentiation of 
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-~ noting that w, but not v and w, is a function of 7, Pr, and 7, 
ey may be put in the following form (with the aid of eq 22): 


oT) _ wvw(vw-- 1) a 
oT, (uvow—1)? (42) 
o7)= _ vw(u+ 1)? 


oT.) ~ (www—1)? (43) 


or, | 
oT ,) (uvew—1 a Mt) 


When eq 42, 43, and 44 are substituted into eq 41, there is obtained 
the relation: 





(2) = (vw 1)?(ute*w? + 1) + ew? (u+ 


(uow—1)* (45) 


which expresses (o;/c) as a function of u, v, and w. Equation 45 
— be put in a more useful form as follows: 

or any given segment of an equilibrium curve, the greatest 
accuracy in 7, (minimum oy) will be obtained by choosing z 2, and 2; 
at the beginning and end of the segment in question (corresponding 
to a maximum value of w). With z, and 2, fixed, the following 
relation then holds: 


uv=c, (46) 


where ¢ is a constant characteristic 7 the given segment of the curve, 
and independent of the choice of z,. Equation 46 may be derived in 
the following way. 

It is evident from eq 22 and 23 that the product (uv) is the ratio 
of the slope of the chord HI to the slope of the chord GH (fig. 2) 
From eq 18 ¢ is then given by 


ii eer) ~ 
e=( Fp , (47) 
which is independent of 2p. 


If one eliminates u from eq 45 and 46, there is obtained the relation 


(2)= _ (vw 1)?(c*w’+ 0?) + w?(e+0)* 
v’(cw—1)* 


A numerical investigation of eq 48 shows, as would be expected, 
that, for the range of conditions encountered in practice, the optimum 
value of » is in the neighborhood of unity (corresponding to z, equl- 
distant between z, and z,). When v=1, w=c, and eq 48 becomes 


(2)- (w+1)*(utw?-++1) +-w*(u+1)* 


(uw—1)* 





(48) 





(49) 


The error in T;,, o;, may be obtained immediately by setting v 
equal to unity in eq 48 or 49. When v=1, 


oy \?_4(u*4-1) + (u+1)* ‘ 
( 1) = Cas) a si 
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From eq 49, (¢;/c), for v=1, has been calculated for w=1, 1.5, 2, 
3,5,and », and for w=1, 1.25, 1.5, 2, 3,4,and ». The results are 
presented in figure 5, in which the abscissas represent (1/w), and the 
ordinates represent (¢;/c), on a logarithmic scale. The curve for 
w=1 represents (cy,/c), and rises to infinity at (1/u)=1, which cor- 
responds to a freezing curve of constant slope. The other curves, 
representing (¢;/0) for various values of w, are finite for all values 
of (1/w). 














fe) 
Figure 5.—Chart to show the uncertainty of the determination of T; and Tso: 


_Figure 5 is expected to be convenient for estimating the uncertain- 
ties of 7’, and 7’, in practice. In calculating T, and 7, from eq 25 


and 33, w and w are first evaluated from eq 22 and 24 (v=1). The 
values of (os/o0) and (,,/7)may then be read immediately from figure 5. 


It remains for the experimenter to estimate a reasonable value for o. 
In any case, figure 5 indicates relative errors for various values of 
wand w. 

The errors discussed above do not include those arising from 
departures from the ideal solution laws or from lack of constancy of 
the rate of crystallization or melting (see appendix, page 212.) 
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VI. ALTERNATIVE METHOD FOR DETERMINING THe 
FREEZING POINT (AND THE FREEZING POINT For 
ZERO IMPURITY) 


The foregoing methods for determining 7; (and, where appropriate. 
T%) may be used when the experimenter is able from previous expe. 
rience to select with some assurance the equilibrium portion of the 
time-temperature freezing or melting curve. When there is some 
doubt as to what portion, if any, of the time-temperature freezing or 
melting curve constitutes the equilibrium portion, the following 
procedure may be used to select that portion of the curve that is 
consistent with the hyperbolic form of eq 16. At the same time, 
there is obtained a value for the freezing pomt, 7;, and within certain 
wide limits, a value for the freezing point for zero impurity, 7h, if 
the observations cover a sufficiently large fraction of the material 
crystallized or melted. 

Let (Tm, 2m) be an arbitrary fixed point chosen on the equilibrium 
portion of a time-temperature curve. Then it may be shown from 
eq 16 that a plot of each of the following functions against the varia- 
ble indicated will yield a straight line: 


(q=7-) versus : ) 
T—T,, 2—Zm 
(= versus (T—T,,), or T 


~—-™ ) versus (2—Z»), OF 2 


~—" 
T— m 


The plot of (2—2m)/(T—T»m) versus (<—Zm) is the most convenient 
for the purposes of this paper. Define new quantities y and z 4s 


follows: 
=(Fo7" =(FoF (51 
y= y = Yr fF 


t=(2—2m); Ly=(Zp—Zm) (52) 


Then it follows from eq 16 that 


Y=Yot sz; Ys=Yot Sry (53 
where 
=( 1 ) (54) 
SAT Tn Pe 
__{_@ Tk’) _ rs 
w=—(r Ta & 


Thus a plot of y against z yields a straight line for that portion 0 
the curve, presumably the equilibrium portion, which is consistetl! 
with eq 16. Extrapolation of the straight line to r=2,(2=%), 
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yields yy; Ty may then be calculated from the equation 


T=Tnt+% (56) 
Yr 
similarly, 
T= Tmt? /8=Tnt 2 ™, (57) 


since tm=0. The smoothed temperature at any other time is given 
by the relation 
PO FO: El (! 
T=T;, ra y 58) 
The point (7, 2m) should be selected near the middle of that por- 
tion of the time-temperature curve which the experimenter regards 
as most likely to represent equilibrium. The linearity of the y versus 
z plot is then an indication of the range on either side of (7, 2m) in 
which equilibrium has apparently been maintained: It should be 
noted that the error in the calculated value of y increases rapidly as 
zdecreases. In plotting y it is better to plot circles, of radii inversely 
proportional to z, than points. This will insure approximately the 
correct weighting of the points in drawing the straight line. It 
should also be noted that a slight displacement of (7, 2m) from the 
true curve will have a negligible effect on the location of the straight 
line, provided (7',, 2m) is about in the middle of the equilibrium por- 
tion of the curve. 
Equation 57 relates 7’, to the slope, s, of the y versus z plot. How- 
ever, it is clear that s is a sort of second derivative of the time-tem- 
perature curve. It is to be expected, therefore, that the value of 7, 


calculated from eq 57 may be.gjuite inaccurate. It is possible to re- 
verse the procedure, if 7, is known, and calculate the value of the 


slope of the y versus z plot to be expected on the basis of eq 16. 

The y versus x plot is principally of use in extrapolating to obtain 
T;, by eq 56, and also in obtaining smoothed values of the tempera- 
ture at any other time, by eq 58. The smoothed values of the tem- 
perature may be used in conjunction with the other methods of 
obtaining 7',, as described [1]. 

(1). 
VII. DISCUSSION 


In connection with the practical application of the foregoing 
methods for determining from time-temperature freezing and melting 
curves the freezing point, 7’, and the freezing point for zero impurity, 
Ty, the following points are important: 


(a) The experimenter must assure himself that the time-tempera- 
ture relations substantially follow eq 16 over that part of the experi- 
ment which represents thermodynamic equilibrium between the 
liquid ae solid phases. (See section VI and the appendix of this 
paper. 

(b) The precision of determining the freezing point, 7, will be 
nearly as good as the determination of temperatures along the equi- 
librium part of the curve, except when the time and temperature 
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covered by undercooling are very large. In terms of the uncertainty 
in the temperatures along the equilibrium part of the curve, the un- 
tainty in 7, is given by eq 49 and illustrated in figure 5. 

(c) The precision of determining the freezing point for zero impw- 
ity, T’;,, by the method described, will be much less precise than the 
determination of the freezing point. In terms of the uncertainty jn 
the temperatures along the equilibrium portion of the curve, the un- 
certainty in 7;, is given by eq 50 and illustrated in figure 5 by the 
curve for which w=1. The practical utility of determining 7, by 
this method depends entirely upon the equilibrium portion of the 
curve covering a sufficiently large range of r (fraction crystallized, or 
melted) so that the product wv in eq 33, the ratio of the slopes of the 
two chords, will have a value appreciably greater than unity. If the 
observations of the given experiment do not cover a sufficiently large 
fraction of material crystallized or melted, no useful determination of 
T;, may be made by this method. For any given experimental pro- 
cedure, it is desirable to check the value of 7;, obtained by this method 
with the value obtained by the method previously described [1], which 
has been shown to give reliable results. 
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IX. APPENDIX 


In the foregoing discussions, it has been assumed for simplicity that the rate of 
crystallization is a constant, or that r is linear in z, as given by eq 15. It is of 
interest to consider the effect of a varying rate of crystallization. Two types of 
variation may be distinguished: 

(a) The rate of crystallization varies initially, but approaches a constant value 
as the experiment proceeds. 

(b) The rate of crystallization varies progressively as the experiment proceeds 
and does not tend to approach a constant value. Supercooling should give rise 
to a variation of type (a). Variations of type (a) will destroy, for a certain time, 
the rectangular hyperbolic relation between 7' and z, eq 16, on which the linearity 
of the (z—z,,)/(T—T,,) versus(z—2z,,) plot depends. It may therefore be assumed 
that in the linear region of this plot such variations have decreased sufficiently to 
be neglected. 

Variations of type (b) are of more importance, since they increase throughout 
the experiment. Such variations would be caused, for example, by any factor, 
such as a decreasing (or increasing) head of temperature, that alters the rate of 
heat transfer and hence the rate of crystallization (or melting). 

If the variation is of type (b), r may be represented by a power series in (z— 2,): 





r=k(z—2z,)[1—e(z—z))+ .. .]. (59) 

Now consider the following expansion, valid for [e(z—z,;)]? less than unity: 
feet |= bee eit 0) 

[; et aaa =k(z—z,)[l—e(z—z-)+ .. .]. (6 


Therefore, when the higher order terms in eq 59 and 60 may be neglected, r may be 
written: 


ae (61) 


r= 


“1+e(z—z) 


oan i GR dee 


a oe me ae, ee ae 
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Equation 61 is a rectangular hyperbolic relation between r and z. Furthermore, 
eq 14 is a rectangular hyperbolic relation between T and r. It is easily shown 
that the relation between 7' and z obtained by substitution of eq 61 in eq 14 is also 
rectangular hyperbolic. Specifically, eq 16 is replaced by the equation: 


i a’e \__a’k'/(k’—e) 
T=(Ty+p-) 1—(k’—e)(z—z))_ (62) 


Equation 62 reduces to eq 16 for e equal to zero. 

t is clear from eq 62 that all the results previously derived on the basis of eq 16 
are unchanged, provided k’ is replaced by (k’—e), a’ is replaced by [a’k’/(k’—e)], 
and Ty, is replaced by [7',+-a’e/(k’—e)]. f the previous results which are unaffec- 
ted, there may be mentioned the theorem concerning the slopes of the chords, eq 
17, and the linearity of the (2—zm)/(7— Tm) versus (2z—Zm) plot. The calculated 
value of 7’ will not be in error. However, the calculated value of Tj, will be in 
error by the amount [a’e/(k’ —e)], corresponding to the fact that the asymptote of 
the hyperbola is shifted. 








WASHINGTON, January 6, 1944. 





U. S. DEPARTMENT OF COMMERCE NATIONAL BurgAu OF STANDARDS 
RESEARCH PAPER RP1586 


Part of Journal of Research of the National Bureau of Standards, Volume 32, 
May 1944 





A METHOD FOR THE DETERMINATION OF THE pH OF 
0.05-MOLAL SOLUTIONS OF ACID POTASSIUM PHTHAL- 
ATE WITH OR WITHOUT POTASSIUM CHLORIDE 


By Walter J. Hamer and S. F. Acree 





ABSTRACT 


The pH values of a 0.05-m solution of acid potassium phthalate containing 
various amounts of potassium chloride were determined at 5-degree intervals from 
0° to 60° C, inclusive, from the measurements of the electromotive force of galvanic 
cells without liquid junction using hydrogen and silver-silver-chloride electrodes. 
A method is described for the determination of the pH directly from the emf by 
means of the equation 


pH=[(£— E°)/(RT/F) +log mci— PJ/Q, 


where P and Q are constants whose numerical values depend on the ionic strength 
of the solution and the nature of the cation and E°, R, T, and F have their usual 
significance. By this method a 0.05-m solution of acid potassium phthalate, 
National Bureau of Standards Standard Sample 84a, is found to have a pH of 
4.008 at 25° C, and this value does not vary much with the temperature. A 
description is also given of the application of the method to the determination of 
the pH and the ionic strength of unknown solutions of low salt content. 

Equations were formulated to express the variation of the pH of a 0.05-m 
solution of acid potassium phthalate with temperature and with the concentration 
of potassium chloride. The pH values may be computed for temperatures from 0° 
to 60° C, inclusive, by the equation 


pH=5.13 log T+ 1519.62/T+0.01092 T — 17.039, 


where T'=t° C+273.16. The pH values for concentrations of potassium chloride 
from 0 to 0.05 m may be computed by the equation 


PH wits salt — PH witho ut salt” 0.993mxc1 + 2. 124mikci. 


The solutions may be readily prepared from known weights of acid potassium 
phthalate, potassium chloride, and distilled water and are well suited for use as 
= standards, in that their pH values do not change much with temperature or 
dilution. 
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I. INTRODUCTION 


Buffer solutions prepared with acid potassium phthalate are ex- 
tensively used for the calibration of commercial pH meters and other 
pH equipment. A 0.05-m solution is generally employed and 

' 215 
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conveniently prepared from the dry salt and distilled water. The 
pH value of this solution has usually been determined from the 
measurements of galvanic cells comprised of various types of hydrogen 
and calomel electrodes. By this method, pH values ranging from 
3.92 to 4.0] [1 to 4] have been reported for a temperature of 25° (. 
The differences in these pH values are due partly to lack of equilibrium 
of the galvanic cell [5, 6, 7], partly to liquid-junction errors [8 to 12 
arising from the use of calomel electrodes, and partly to the use of 
hydrogen-ion concentration rather than hydrogen-ion activity as the 
basis of the definition of the pH scale.? This paper gives a method for 
the determination of the pH of a 0.05—m solution of acid potassium 
ey from the measurements of the,emf of galvanic cells without 
iquid junction. The effect of varying amounts of potassium chloride 
upon the pH of a 0.05—m solution of acid potassium phthalate is also 


.given. The latter data are included because they give information 


concerning the activity coefficients of the ions composing the mixture: 
A description is also given of the application of the method to the 
determination of the pH values and ionic strengths of unknown solu- 
tions of low salt content [24]. 


II. EXPERIMENTAL PROCEDURE 


The method is based upon the measurements of the emf of the 
galvanic cell 


Pd, H2(g)|KHPh(0.05 m), KCl(m)|AgCl(s)|Ag(s) 


for a series of concentrations (m) of potassium chloride, in which 
KHPh designates acid potassium phthalate, and the other symbols 
have their usual significance. 

The hydrogen electrodes used were made of platinum foil (10 by 
15 by 0.2 mm), and were plated with a thin coat of palladium sponge 
by electrolysis of a 1-percent solution of palladium chloride contain- 
ing a trace of lead acetate for 1 minute at a current density of 0.04 
amp/cm? [5]. Electrodes of this type do not catalyze the reduction 
of acid potassium phthalate and function reversibly for well over 100 
hours in solutions containing various amounts of potassium chloride. 
Silver-silver-chloride electrodes of the thermal-electrolytic type were 
prepared from pure materials, as described in a previous paper [13]. 
They were sufficiently aged to bring them to a stable state. Other 
equipment, including the cells, thermostat, and emf-recording instru- 
ments have been described in detail [13]. 

The acid potassium phthalate used was National Bureau of Stand- 
ards Standard Sample 84a, having a certified purity of 100.00 percent, 

1 sy ape in brackets indicate the literature references at the end of this paper. 

2 It is well known that an exact differentiation between hydrogen-ion concentration and activity cannot 
be made, since the value of the activity coefficient of the hydrogen-ion cannot be determined without some 
assumption. Hitchcock and Taylor (3, and MacInnes, Belcher, and Shedlovsky [4] in their methods obtain 
8 quantity which approximates the value of the hydrogen-ion activity. In the work on pH standards in this 
Bureau the assumption is made that the activity coefficient of the hydrogen-ion is equal to the mean activity 
coefficient of hydrochloric acid, a quantity which can be determined experimentally from electromotive: 
force a See references {8, 13 to 15] for discussions of this assumption in connection with work 
on 8 ards. 

if the dissociation of an acid salt conformed to the simple law of mass action for finite concentrations, 
the effect of a neutral salt on the dissociation of the acid salt could be determined solely from its concentra 
tion. However, asa result of ionic association, the nature of the medium, and salt effects, en tes 
occur which lead to departures from ideality. For a review of the earlier studies of the effects of salts on Lad 


dissociation and reactions of weak electrolytes, see the paper by Burton, Hamer, and Acree [16]. See 
the papers by Loomis and Acree [17] and Bronsted [18]. 
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determined by titration. It was prepared for use according to the 
directions given in the Bureau certificate. Potassium chloride was 
recrystallized twice from conductivity water and then thrice by the 
addition of 95-percent ethanol to remove the last traces of bromide. 
It was dried for at least 10 hours at 110° C and allowed to cool in a 
desiccator, using concentrated sulfuric acid as the desiccant. The 


2] solutions were prepared from known weights of the dry salts and 
of conductivity water. The concentrations were expressed in terms of 
le weight of solute per 1,000 g of water. All the solutions were deaerated 
or by means of purified hydrogen gas before they were introduced into 
n the cells. Corrections were made for the loss of water from the solu- 


tions during the deaeration process and for the differences between the 
densities of hydrogen and air. The specific conductance of the water 
used in this work was generally 1.0X10~® mho and sometimes was as 
low as 0.6 X10~° mho. 

In table 1, the values of emf given were corrected in the usual 
manner to 1 atmosphere of hydrogen pressure. The results of tripli- 
cate (three) determinations (values) are given for 25° C and of dupli- 
cate determinations (two) for 45° C. The values of the emf were first 
recorded at 25° C. The bath was then cooled to 0° C in 5-degree 
steps and the emf recorded at each temperature. Then the tempera- 
ture of the bath was raised to 25° C and the emf again recorded. The 
temperature of the bath was then raised to 60° C in 5-degree steps 
and the emf recorded at each temperature. Finally, the bath was 
cooled first to 45° C and then to 25° C and the emf measured at each 
temperature. The first and final emf at 25° C agreed within 0.09 mv, 
or 0.002 pH unit, on the average for the nine concentrations. 


TaBLE 1.—Electromotive forces of the galvanic cell Pd, H2(g)|KHPh(0.05 m), 
KCl (m)|AgCl (s)|Ag (s) from 0° to 60° C 


























y 
e {Hydrogen electrodes were aged for 1 hour; silver-silver-chloride electrodes were aged for 12 hours.] 
Molality of potassium chloride=m 
) Tempera- 
ture 
0.05 0.03 0.02 0.01 0.007 0.005 0.003 0.002 0.001 
°C v 7] D v v v v v v 
aveceest Qe 0. 58976 0. 54962 0. 56611 0. 57459 0. 58255 0. 59476 0. 60420 0. 62059 
sided fh te . 52976 . 54255 . 55261 . 56936 . 57803 - 58620 . 59849 . 60812 . 62478 
Leckie nan . 53229 . 54532 . 55550 . 57263 . 58144 . 58963 . 60222 . 61210 . 62908 
ixtiocwend . 53482 . 54806 . 55847 . 57583 . 58485 . 59322 . 60591 . 61602 . 63334 
SS . 55085 . 56143 . 57915 . 58828 . 59678 . 60968 . 61994 . 63760 
We kia cia . 53999 - 55363 . 56440 . 58255 . 59173 . 60040 . 61354 . 62396 . 64186 
, ge . 53998 . 55364 - 56435 . 58257 . 59174 . 60042 . 61356 . 62393 . 64190 
We cone - 54003 . 55369 . 56438 . 58270 . 59184 . 60045 . 61368 . 62403 . 64201 
Cae, ee, . 55644 - 56738 . 58566 . 69514 . 60394 . 61731 . 62789 . 64618 
. 54510 . 55919 - 57023 . 58888 . 59854 . 60748 . 62107 . 63187 . 65036 
eS . 54752 . 56190 . 57318 . 59206 . 60187 . 61096 . 62475 . 63571 . 65457 
GBF iscsicinee . 55008 . 56452 . 57600 . 59528 . 60514 . 61443 . 62847 . 63963 . 65866 
ke . 55011 . 56455 . 57606 . 59529 . 60519 . 61454 . 62853 . 63972 . 65875 
tee, . 55264 . 56746 . 57908 . 59859 . 60871 . 61810 . 63234 . 64366 . 66304 
poeesesce . 55519 . 57023 . 58207 . 60184 - 61210 . 62164 . 63610 . 64767 . 66729 
pene cunts . 55773 . 57287 . 58487 . 60506 . 61543 . 62508 . 63974 . 65152 . 67137 




















*emf recorded 6 hours after the start of flow of the hydrogen gas. 

>emf recorded after measurements at 0°, 5°, 10°, 15°, 20°, and 25° (first) C. 
‘emf recorded after measurements at all the other temperatures were made. 
‘emf recorded after measurements at the lower temperatures were made. 
*emf recorded after measurements were made at 50°, 55°, and 60° C. 
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III. METHOD FOR THE DETERMINATION OF THE pH OF 
ACID POTASSIUM PHTHALATE 


The relation between the emf of the cell and the hydrogen- and 
chloride-ion activities of the solution is given by the equation 


E= E°—(RT/F) In (dg¢q1) = E°— (RT/F) In (fafarmamar), (1) 


where E is the measured emf; E° the potential of the silver—silver- 
chloride electrode at unit activity of hydrochloric acid; a, f, and m 
are, respectively, the activity, activity coefficient, and molality of the 
ionic species denoted by subscripts; and R, 7, and F have their usual 
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Figure 1.—Plots of the numerical values of —log (fufcimn) at 0° to 60° C for solu- 
— af 0.05-m acid potassium phthalate containing different amounts of potassium 
chloride. 


0.09, 0.06, 0.04, 0.03, 0.02, and 0.01 were subtracted from the values at 0°, 5°, 10°, 15°, 20°, and 25° C for 
convenience in plotting. 
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significance. Upon collection of all known quantities on the left, eq 


1 becomes 
(E— E°)/(RT/F) +log ma=—log (fafcimna), (2) 


by which values of —log (fafcimg), a quantity proportional to pH or 
—log (fama), may be calculated. It would therefore be a simple 
matter to calculate pH values from the emf if values of —log fq for 
the solution were known. Hamer and Acree [13] have shown how 
values of —log fo, may be determined from the emf of cells without 
























” liquid junction, using the equations of Debye and Hiickel. This 
: method entails comprehensive studies, however, of the thermody- 
a 


namics of the system under study, in this case of various mixtures of 
o-phthalic acid and potassium hydroxide. However, a simpler 
method, although probably less precise, may be followed and is herein 
described for the temperature 25° C. 

From the emf determinations recorded in table 1, values of—log 
(fafciMm) Were calculated by eq 2 for each concentration of potassium 
chloride for the different temperatures, using the data for E° deter- 
mined by Harned and Ehlers [19] and values of R, 7, and F employed 
by Hamer, Burton, and Acree [14], and are given in table 2. These 
values are then plotted against the molality of potassium chloride, 
as shown in figure 1, and extrapolation made to zero salt concentra- 
tion, Which gives the value of —log Un foima)” for 0.05-m acid potas- 
sium phthalate in the absence,of chloride ion. The superscript, zero, 
is used to denote quantities for solutions containing no chloride ion. 
Within the experimental error (0.002 in pH unit), the values vary 
linearly with the molality of potassium chloride. This variation may 
be expressed by 


—log (fafoima)= —log (fafcima)? —kmxc, (3) 


where k is a constant whose magnitude is given by the slope of the 
straight lines of figure 1. Values of —log (fafcymx)° and k for the 
different temperatures were determined by the method of least 
squares and are given in the last two columns of table 2. At 25° C, 
—log (fafoimn)° equals 4.0934. 


TABLE 2.—Numerical values of —log (fuafcimn) in 0.05-m acid potassium cs 
containing various amounts of potassium chloride from O° to 60° 





























Molality of potassium chloride 

Tempera- 

ture ' | 

0.05 | 0. 03 | 0. 02 | 0.01 | 0.007 | 0. 005 | 0.003 | 0.002; 0.001 | 0. 000 | kb 

°C | 
0............| 4.0665 | 4.0755 | 4.0814 | 4.0847 | 4.0865 | 4.0871 | 4.0906 | 4.0888 | 4.0903 | 4.0903 | 0.48 
Cars tacit 4.0706 | 4.0768 | 4.0794 | 4.0816 | 4.0836 | 4.0845 | 4.0830 | 4.0840 | 4.0851 . 48 
10 | 4.0581 | 4.0683 | 4.0735 | 4.0775 | 4.0794 | 4.0792 | 4.0815 | 4.0813 | 4.0826 | 4.0826 .49 
5 fe 4.0596 | 4.0695 | 4.0756 | 4.0782 | 4.0811 | 4.0815 | 4.0817 | 4.0825 | 4.0845 | 4.0840 .48 

PREIS ie: 4 4.0731 | 4.0790 | 4.0827 | 4.0848 | 4.0850 | 4.0850 | 4.0855 | 4.0880 | 4.0875 48 
oS 4.0695 | 4.0783 | 4.0843 | 4.0902 | 4.0904 | 4.0909 | 4.0913 | 4.0915 | 4.0930 | 4.0934 48 
Se 4.0772-| 4.0872 | 4.0930 | 4.0960 | 4.0987 | 4.0990 | 4.0004 | 4.0995 | 4.1025 | 4.1016 48 
Re 4.0980 | 4.1026 | 4.1067 | 4.108 | 4.1100 | 4.1105 | 4.1110 | 4.1125 | 4.1122 46 
ne ee 4 4.1095 | 4.1151 | 4.1180 | 4.1210 | 4.1213 | 4.1215 | 4.1219 | 4.1245 | 4. 1237 47 
458 4.1167 | 4.1238 | 4.1291 | 4.1340 | 4.1355 | 4.1866 | 4.1372 | 4.1379 | 4.1385 | 4.1386 45 
WU iicwinewwacea 4.1338 | 4.1432 | 4.1485 | 4.1518 | 4.1548 | 4.1552 | 4.1555 | 4.1560 | 4.1575 | 4.1574 47 
RE ES 4.1508 | 4.1600 | 4.1659 | 4.1685 | 4.1712 | 4.1717 | 4.1720 | 4.1736 | 4.1740 | 4.1741 - 46 
60...........] 4.1702 | 4.1775 | 4.1830 | 4.1875 | 4.1895 | 4.1895 | 4.1895 | 4.1916 | 4.1910 | 4.1917 44 





























* Only the first recorded emf were used at these temperatures. Data obtained with the second or third 
measurements of the emf agree with the above data within0.00090n the average. kis the coefficient in the 
equation —log (fafcrmn) =—log (fafcima)°— kmxci where —log (fafc:mn)° represents the quantity —log 
(fafcimg) for zero salt concentration. 








220 Journal of Research of the National Bureau of Standards 


The pH value of the solution may then be obtained from this 
value by use of a large-scale plot of —log (famx)° against —log 
(fafcimn)° for the same ionic strength, or determined from a more 
sensitive plot of the differences between —log (famg)° and —log 
(fufcimn)° against —log (fafcima)°. <A description of the method 
used in the construction of the latter plot follows. Similar measure. 
ments are made with 0.05-m hydrochloric acid and 0.05-m potassium 
hydroxide. This alkali is used, since acid potassium phthalate js 
considered here; for sodium or lithium salts, sodium or lithium 
hydroxide would be used. From the date of Harned and Ehlers 
[19, 20] and of Harned and Hamer [21], values of 1.4628 and 12.7021 
are obtained for —log (fx fcims)° for hydrochloric acid and potassium 
hydroxide, respectively. Also by means of the equations 


—log(/amy) =pH (4) 


—log Kyw/(fou™on) = —log (fumy) =pH, (5) 


the pH values of pure hydrochloric acid (eq 4) and pure potassium 
hydroxide (eq 5) may be calculated, 3° lay is the ionization con- 
stant of water and m and f are the molalities and activity coefficients 
of the ions denoted by subscripts. Values of 1.3817 and 12.6114 are 


and 
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Figure 2.—Plots of the differences between —(log fumu)° and —log (fafcimn)’ 
against —(log fafcimu)° at 25° C for various ionic strengths of hydrochloric acid 
and potassium hydrowide. 


Data for acid potassium phthalate are read from = plots, as shown by the dotted vertical and hori 
zontal lines. 


AALN Dt ee a ~ , ae, | ay =.4°e = i tes bs ae a. 
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obtained by these equations for the pH of 0.05-m solutions of hydro- 
chloric acid and potassium hydroxide, respectively, using 0.8304 for 
fa [19], 1.008107" for Kw [21], and 0.824 for fon [22, 20]. Hence 
—log (fufcimu)°® and —log (famx)° differ by 0.0811 for 0.05-m hydro- 
chloric acid and by 0.0907 for 0.05-m potassium hydroxide at 25° C. 
If the assumption is made that a straight line may be drawn between 
the points located by plotting these differences against —log (fu fovma)° 
for a definite ionic strength, it is then possible to determine the 
difference between —log (fufcima)° and —log (famg)° for 0.05—m 
acid potassium phthalate from the plot.® 

In figure 2 these differences are shown plotted against values of 
—log (fafcima)° for ionic strengths of 0.02, 0.05, 0.10, and 0.20 
calculated in a similar manner from data in the literature [19 to 22] 
for 25° C. From the plot, a difference of 0.0831 is read for 0.05—m 
acid potassium phthalate, for which —log (fyfcymg)° is 4.0934; hence 
its pH value is 4.0934—0.0831, or 4.010 at 25° C. 

However, the ionic strength of 0.05-m acid potassium phthalate is 
somewhat larger than 0.05, owing to several ionic reactions which 
occur in an aqueous solution of this salt. When acid potassium 
phthalate is dissolved in water it ionizes according to the equation 


KHPh->Kt+HPh-, (6) 
and since HPh~ behaves as an acid, it ionizes further: 
HPh-=Ht+ Ph", (7) 


to give a certain amount of the bivalent phthalate ion. However, 
H* does not equal Ph*, since part of the H* ions reacts with HPh~ 
according to the reaction 


H++HPh-=H,Ph, (8) 


to give a definite amount of undissociated o-phthalic acid, H,Ph. 
The amounts of H,Ph, HPh~, and Ph™ may be calculated by using 
the mass action equations for the above reactions and for the two 
steps in the ionization of o-pbthalic acid with the aid of the Debye- 
Hiickel equations for the activity coefficients. These were found to 
be, respectively, 0.0032, 0.0436, and 0.0032, using 0.0013 and 0.0000038, 
respectively, [23] for A, and Ko, the first and second ionization con- 
stants of o-phthalic acid. Therefore, the ionic strength of a 0.05-m 
solution of acid potassium phthalate is 0.0533 instead ® of 0.05. In 
figure 2, the dotted line corresponds to this ionic strength. The posi- 
sition of the dotted line was determined from that of the others by 
means of plots of [log (famg)°—log (fafcimn)°] for pH values of 0 
and 14 against the ionic strength. A difference of 0.0854 is then read 
from this line for acid potassium phthalate, and hence its pH value 
is 4.0934—0.0854, or 4.008, which is only 0.002 pH unit lower than 
the value based upon the stoichiometric ionic strength end is the 
value certified by this Bureau for ecid potassium phthalate Standard 
Sample 84a. 

‘The assumption is made that for for=-+/x/on for the alkaline solutions and fa=~fafc for the acid 
solutions, or that the activity coefficient of an ion in an electrolyte is equal to the mean activity coefficient 
of the ions composing the electrolyte. This does not imply that fu=fx=fci=fon. See E. A. Guggenheim, 
J. Phys. Chem. 84, 1758 (1930) for a discussion of this assumption. 

' The small contribution of mx to the ionic strength can be calculated from the provisional values of pH. 

§ The ionic strength is equal to 1/2 Zimjz;*, where mj and 2; represent the molality and valence of each 


A A vies. For the above case the ionic strength equals mypn+3mpn+1/2mH or (0.0436-+-0.0096-+0.00005 = 
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This method is estimated to be accurate within 0.005 to 0.01 pq 
unit, even if the stoichiometric ionic strength is employed instead of 
the true ionic strength.’ The pH values based upon en ionic strength 
of 0.05 and of 0.0533 are given in columns 3 and 4 of table 3 for the 
temperatures * 0° to 60° C. Data obtained by former investigator 
by other methods [1 to 4] are given for comparison. The values ob. 
tained in this investigation are 0.04 to 0.05 unit higher than those of 
Clark and Lubs [1] and Kolthoff and Tekelenburg [2], who used the 
concentration instead of the activity of hydrogen ions as a basis for 
their pH scale, but are in close agreement with the values of Hitch- 
cock and Taylor [3] and of MacInnes, Belcher, and Shedlovsky {4}, 


TaBLeE 3.—pH values of 0.05—m acid potassium phthalate at temperatures from 
0° to 60° C, inclusive 


{Data of other investigators are given for comparison] 











Clark and} Kolthoff and | Hitchcock | MacInnes, 
‘ an Arner Jan ged Lubs | Tekelenburg | and Taylor | Belcher, and 
Temperature) —log(fufcimn) strenght of | strength of (1916) (1927) (1937-38) ~~ 
°C 
0 4. 0903 4. 009 ft SESS SERA Se pape ee ie 
5 4. 0851 4.004 EN os nacdwndwledas caneguctadalaestGeabmuneticeiiiettniac: 
10 4. 0826 4.001 OS See See ye em peep epee. 
12.5 | Be SEAN NRE SN eae PEE Sibi ottien nc ewe knancgenSUlbwucwesteuaee 4. 000 
15 4. 0840 4. 002 Dongs erebeinc se penandceedlwndnccemebrlibipesiikenncs 
ee See er ee Aopen ©4001 f..~42555- ©3048 CO 4....-5 ac 
20 4. 0875 4. 005 4, 003 SF RNS, CERT ARS Oe OER a 
25 4. 0934 4.010 hd AE? St RARE IS ELE 4. 008 4.000 
30 4. 1016 4.018 rE od Cee BGG ASOD | bois Satcdcwedsbobenccic. 
35 4. 1122 4. 027 SRR ener: Men Eee. a an 
WS Btiwenncanton i eer Ale epee CEE Be oe cdpancctens ehh ations 4. 025 4.015 
40 4. 1237 4. 038 O.00E 1s d+ on 5 Xe |e pee Rae ie 
45 4. 1386 4.052 Ng be SER RRR SESS ee Beare AB ED SES ae 
50 4. 1574 4.070 St eee ek eee ee See ee 
55 4. 1741 4. 086 DE Ac cbedn oh giledsdcadeoudddl bie cnecs 
60 4.1917 4. 102 fi eee 2S 4.05 (6.06 fo. 5.5 2ul5e. oO eS te 





























— ®.. 
* pH values were calculated by the equation pH= (og fale ma £ + 

> Values calculated by equation 13. 

® Values in parentheses were determined with a beep a electrode, whereas the other values of Kolthoff 
and Tekelenburg were determined with a quinhydrone electrode. They considered the hydrogen-electrode 
values less reliable because of possible reduction of the acid potassium phthalate. 


The straight lines of figure 2 may be expressed by the equation 


log (fumuz)°—log (fufcimu)°=a—b log (fafoimn)®, (9) 


where a is the value of log (famg)°—log (fafcimn)® for —log 
(fafcima)°=0, and 6 is a constant whose magnitude at various ionic 
strengths is given by the slopes of the straight lines. This equation 
may be converted to the form 


pH=[—log (fafcimu)°— P)/Q, (10) 


7 It is difficult at the present time to give an evaluation of the accuracy of this method. Further work 
on the thermodynamics of mixtures of o-phthalic acid and potassium hydroxide (publication pending) will 
shed light on this subject. However, the close agreement of the pH value for 0.05-m acid potassium 

hthalate at 25° C obtained by this method with the values obtained by Hitchcock and Taylor and by 

acInnes, Belcher, and Shedlovsky given in table 3 indicates that this method is probably accurate 
within 0.005 to 0.01 pH unit. ‘ 

8 Harned and Cook [11] obtained activity coefficients at only 0°, 10°, 20°, 25°, and 35° C, Values at the 
other temperatures were obtained by extrapolation. 
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whire P and Q are new constants. Values of P and Q for ionic 
strengths of 0.05 and 0.0533 are given in table 4 for temperatures from 
0° to 60° C, inclusive. Values of P and Q for different ionic strengths 
at 25° C are given in table 5. Furthermore, since —log (fafcimg) 
equals (E—E*°)/(RT/F)+log meu, eq 10 may be written 


pH=[(Z— E°)/(RT/F) +log ma— P)/Q, (11) 


by which the pH may be calculated directly from the emf and the 
numerical values.of P and Q given in table 4, remembering that 
P and Q depend upon the ionic strength and are for solutions contain- 
ing no chloride ion. For solutions containing Cl~ ions, different 
values of P and Q must be used, which are given in table 6. 


TaBLE 4.—Numerical values of the constants P and Q of equation 10 for ionic 
strengths of 0.05 and 0.0533 














Ionic strength 
Temper- ee 
ature | 
0.05 0.0533 

“ P Q od Q 

0 0.0772 1.00107 0.0787 1.00112 

5 . 0773 1.00102 . 0789 1.00106 
10 . 0776 1.00096 - 0795 1.00101 
15 . 0782 1. 00094 . 0799 1. 00098 
20 . 0789 1.00090 . 0805 1.00094 
25 . 0796 1.00085 . 0814 1.00089 
30 . 0807 1. 00082 . 0824 1. 00086 
35 . 0818 1.00079 . 0835 1. 00083 
40 . 0829 1.00076 . 0847 1.00080 
45 . 0838 1. 00073 . 0856 1.00076 
50 . 0849 1.00068 . 0867 1.00072 
55 . 0859 1.00066 . 0877 1.00070 
60 . 0869 1.00065 . 0888 1. 00067 























TaBLE 5.—Numerical values of the constants P and Q of equation 10 for different 
ionic strengths at 25° C 











Ionic 
strength P Q 

0. 005 0.0313 1. 00040 
01 . 0423 1.00055 
02 . 0575 1. 00068 
05 . 0796 1. 00085 
10 . 0976 1. 00092 
20 . 1146 1.00096 

















The pH values of the solutions of 0.05-m acid potassium phthalate 
containing various amounts of potassium chloride were calculated 
similarly, using the proper values of P and Q. The values are given 
in table 7, and in figure 3 are shown plotted against the molality of 
potassium chloride for temperatures of 0°, 25°, 45°, and 60° C. The 
plots exhibit a slight curvature, appearing at molalities between 0.03 
and 0.04. The variation of pH with salt concentration may be ex- 
pressed at all the temperatures by the equation 


PH witn salt — PH witnout sat — 0.993Mx01+2.124m'* Kor. (12) 
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In figure 4 the pH of 0.05-m acid potassium phthalate (curve A) 
and of 0.05-m acid potassium phthalate containing 0.02-m potassium 
chloride (curve B) and 0.05-m potassium chloride (curve C) are shown 
plotted against the temperature in degrees centigrade. The curves 


TaBLE 6.—Numerical values of constants P and © for various concentrations of 
potassium chloride in 0.05-m acid potassium phthalate for use in equation 10 


= 








Molality of potassium chloride 


























Tempera-| 0.001 | 0.002 | 0. 003 | 0.005 | 0.007 0.91 0.02 | 0.03 | 0.05 
P VALUES 
l 
°C | 

| 
0} 0.0783} 0.0788| 0.0794} 0.0805 | 0.0811} 0.0830] 0,0872| 0.0016 | 0.09% 
5| .0701| .0706| .0804|  .0814/ .0820| .0838|  .0880|  .0926 “0088 
10} :0709}  :0804} .0810| .0822|  .0829|  .0847| 0800}  . 0036 - 0999 
15| 10807} 0813} .0819|  .0820)  .0837|  .0856 | 0809) ©. 0046 ‘1000 
20| 0810} 0821} —. 0827 | +0889 0845 | 0864} 0909 0956 1010 
25 | .0824| .0830|  .0836| 0848 0855 0874} 0919 0965 1029 
30; 0833 0836; .0841|  .0857 0865 0928 0975 1042 
35| .0842| 0847} 0853 | 0866 0874 0892 | 0938 0985 1052 
40| [0851 | :0857/ .0864| .0875| .0883|  .0903|  :0948| [0996 ‘1065 
45} .0861| 0867) .0873| .0885| .0804| 10913} 10950) 1008 ‘1076 
50|  . 0871 .0877 | .9884 | .0806| 0904) .0024| 0971} .1020 . 10% 
55| 0881 10887} .0804/ 0006! :0015| 10035; :0981/ |1032 "1101 
60 | .0801, .0808, .0904' .0917 .0923' 0046)  .0904| | 1044 Wid 








1.00108 | 








0 1.00111 | 1.00110 | 1.00110 1, 00109 | 1, 00109 | 1.00103 | 1.00097 | 1, OOORY 

5 1.00105 1.00104 | 1.00104 1.00103 | 1.00103 | 1.00102; 1.00097 | 1.00091 | 1, 00082 
10 1.00100 1.00099 | 1.00099 | 1.00097 | 1.00097 | 1.00096 1.00091 | 1.00084 1, 00075 
15 1. 00097 1. 00096 1.00096 | 1.00094 { 1.00093 | 1, 00091 1, 00085 1. 00078 1, 00068 
20 1. 00093 1. 00092 1. 00092 | 1.00090 | 1.00089 | 1.00087 1. 00081 1, 00073 1. 00063 

| 
25 1. 00088 1. 00087 1. 00087 1.00084 | 1.00084 | 1. 00082 1. 00076 1. 00068 1, 00058 
30} 1.00085 | 1.00084 1.00084 1.00082, 1.00081 | 1, 00079 1. 00073 1. 00065 1, 00055 
35 1.00082 | 1.00081 1.00081 | 1.00079 | 1.00078 | , 1.00076 1. 00070 1. 00062 1, 00052 
40 1, 00079 1. 00078 1.00078 | 1.00076 | 1.00075 | 1.00073 1, 00067 1. 00059 1, 00045 
45 1, 00075 1, 00074 1. 00073 1, 00071 1, 00069 1, 00067 1. 00061 1, 00053 1, 00048 
50 | 1.00071 1.00069 | 1.00068 | 1.00066 | 1.00064 | 1.00062 | 1.00055 1. 00046 1. 00035 
55 1. 00069 1. 00067 1.00066 | 1.00064 | 1.00062 1. 00060 1. 00053 1. 00044 1. 00033 
60 | 1. 00066 1, 00064 1.00063 | 1.00061 | 1.00059 1, 00057 1. 00050 1. 00041 1, 00030 
} | i i 

















TABLE 7.—pH values of 0.05-m acid potassium phthalate containing various amounts 
of potassium chloride at temperatures from O° to 60° C, inclusive 





Molality of potassium chloride 























Tempera- ier eae 

ture | j 

0.001 | 0.002 | 0.003 | 0.005 | 0.007 | 0.01 | 002 | 0.03 | 0.05 
| 

"e rs be a 
po tiiscies 4.008 | 4.006| 4.007} 4.002] 4001] 3.997) 3.990} 3.980/ 3.965 
ates ee 3.900 | 4.000} 3.998] 3.996] 3.992| 3.985| 3.974 3. 950 
10} «= 999 | «= 3.997} = 3.907] = 3.903 | «3.993 | 3.980) 3.981! 3.971! 3.955 
15........} 4.000] 3.997] 3.996] 3.905) 3.994] 3989) 3.982, 3.972 3, 956 
20..--"-27:] 4.008 | 4.000) 3.999} 3.908} 3.997) 3.998} 3.985) 3.975 3, 960 
~ Rae 4.007| 4.005| 4.004/ 4003; 4002] 4.000; 3.989, 3.970! 3.964 
30...) «016 | «= 412 | 4.012 | = 010 |) 4.009 | += 4.005 | += 3.997} 3.987} 3.971 
35. “| 4025 | 4023 | 4.022| 4.020, 4019] 4014] 4.006| 3.907) 3.982 
40. “| 4036} 4.033} 4032) 4.031! 4.030) 4025|/ 4018) 4.008| 3.902 
45 “| 4049} 64.048) 394.047 4.045) 394.043] 94.040] 4.031 | 4021) 4.008 
50 4.068! 4.066) 4.064, 4063; 4062] 4.057! 4049/ 4.039; 4.02 
55 4.083! 4.082 4.083 4.079| 4.077; 4.073 4.066! 4.055 4.040 
60 4.099 4.009, 4.097) 4.095 ron 
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are displaced a constant amount at all temperatures by the addition 
of potassium chloride and exhibit » minimum at about 10° C. The 
variation of the pH of a 0.05-m solution of acid potassium phthalate 
containing no potassium chloride may be represented by the equation 


pH=5.13 log T-+1519.62/7'+-0.01092 T— 17.039 (13) 


where 7=t°C+-273.16. For the solutions containing salt, the equa- 
tion is modified by changing 17.039 by an amount determined by 
eq 12. The pH values are constant within 0.005 pH unit between 
(° and 30° C, and this fact increases their usefulness as pH standards. 

The potentiometric method herein described is a general one, in 
that it may be used equally well for other concentrations and other 
buffer solutions for the temperature range of 0° to 60° C. It has the 
advantage of being less involved, although less precise, than ones 
based upon a complete thermodynamic study of the system, and the 
pH value of a solution may be determined directly from the electro- 
motive force of galvanic cells without liquid junction. 
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FicurE 3.—Variation in the pH of 0.05-m acid potassium phthalate with the molality 
of added potassium chloride at 0°, 25°, 45°, and 60° C. 
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The method may also be used to determine approximately the pH 
values and ionic strengths of unknowns [24] of low salt content, but if 
chloride salts are present, the chloride concentration must be determined 
for use in the equations. The computations are simpler if the un- 
known contains no chloride. - First, two different amounts of potas- 
sium chloride are added to two separate portions of the unknown, and 
—log (fa fcimu) is calculated by eq 2 for each portion from the elec- 
tromotive force and the known chloride-ion concentration. The 
value of —log (fafcima)° for the solution without added chloride is 
then obtained by means of a straight-line plot of the values of —log 
(fafcimx) against the chloride-ion concentration. This value js 
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FicurE 4.—Variation in the pH of aqueous solutions of mixtures of acid potassium 


phthalate and potassium chloride. 


A, 0.05-m acid potassium phthalate; B, 0.05-m acid potassium phthalate +0.02-m potassium chloride: 
C, 0.05-m acid potassium phthalate +0.05-m potassium chloride. 





then converted into pH units by means of figure 2 or by eq 10. How- 
ever, the ionic strength of the unknown must be determined in order 
to decide which line of figure 2 or what values of P and Q of eq 10 to use 
to obtain the pH value. The ionic strength of the unknown may be 
obtained from the two measured values of —log (fa fcmg) by the 
equation 


—log (fa fewmn):t+log (fafcwms)2= 





1.0196 Vz-+e, 1.0196 Vz-++e, 
1+ yz+a 1+ yz+e, 


where the subscripts 1 and 2 refer to the two different chloride-ion 
concentrations (c; and c) that were employed and z is the ionic strength 
of the unknown. This equation is applicable for ionic strengths of 0.2 
or less. An illustration of the method follows. 

An unknown solution ® upon the addition of 0.01—-m and 0.02-m 


(14) 


* This solution contained 0.01 mole of acid potassium phthalate and 0.01 mole of dipotassium phthalate 
per 1,000 g of water. 
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potassium chloride gave, respectively, electromotive forces of 0.65431 v 
and 0.63472 v. These give, respectively, 5.3035 and 5.2743 for —log 
(fafoimu). The difference (0.02922) of the two values added to 
the first value gives 5.3327 for —log (fafcimu)° for the unknown solu- 
tion containing no chloride ion. The same difference substituted in 
eq 14 gives 0.0037 for x, the ionic strength of the unknown. From 
interpolation of the data of table 4, P and Q are, respectively, 0.0200 
and 1.00025. Therefore, by means of eq 10 the pH of the unknown 
is (5.3327 —0.0200)/1.00025=5.311. This solution was known to 
have an ionic strength of 0.004 and a pH of 5.326, which agree ade- 
quately with the calculated values. ; 

Similar applications to other unknowns may not be as successful. 
However, the method may be employed with a precision of about 
0.1 or better in pH for solutions of ionic strength less than 0.2. The 
values of P given in table 5 correspond to the errors in pH at different 
ionic strengths if no corrections were made for the ionic strength. 
For example, in the above illustration the error in the pH of the 
unknown solution would be 0.02 if no corrections had been made for 
the ionic strength. Furthermore, the method cannot be applied 
successfully to a chloride solution unless its concentration is known 
because the emf is a function of the activity of the chloride ion. 
Additional studies of effects of polar and nonpolar salts of high con- 
centrations on the pH of solutions should make possible an extension 
of the method and theory to a wider range of ionic strengths. 
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THERMAL PROPERTIES OF MOIST FABRICS 
By Charles W. Hock, Arnold M. Sookne, and Milton Harris ! 


ABSTRACT 


The “chilling effect’’, or “clamminess’’, that moist fabrics produce when in 
contact with the body was evaluated by subjective tests, by measurement of the 
drop in temperature that ensued when the moist fabrics were placed on an arti- 
ficial “skin”? surface, and by tests with a moisture-sensitive paper designed to 
measure the extent of contact which the fabrics made with a surface. Using 
fabrics of various fiber compositions and constructions, a good qualitative rela- 
tion was found in these tests. Fabrics which produced considerable chilling in 
subjective tests were found to make good contact and to cause a substantial drop 
in skin temperature. Conversely, fabrics which caused little or no clamminess 
made poor contact and the accompanying drop in temperature was relatively 
small, The results of these experiments show clearly the progressive improve- 
ment of the fabrics with respect to chilling, as their wool content is increased, and 
also the superiority of certain types of construction that minimize the extent of 
contact of the fabrics with the skin, 
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I. INTRODUCTION 


Practical experience has demonstrated that moist fabrics in contact 
with the body produce an unpleasant sensation, commonly referred 
to as a “chilling effect’? or a “clammy feel.” The intensity of the 
sensation varies with different fibers and fabrics; thus, on the basis of 
general experience, the merits of wool fabrics over similar cotton 
fabrics have long been recognized. This is one of the reasons, for 
example, for the preference usually given to woolen underwear and 
other garments for use in cold climates under conditions where 
physical labor causes considerable perspiration. It also accounts for 
the preference commonly expressed for woolen bathing suits. 

In selecting fabrics for warmth, the thermal resistance is a factor 
of primary importance, and this property has been measured by 
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many investigators [1, 2,3, 4, 5, 6].2_ In general, the results indicate 
that for comparable thicknesses and densities there are no outstanding 
differences between different fabrics. The chilling effect noted above, 
however, has suggested that marked differences may be obtained as 
the amount of moisture contained in a fabric is increased, and 
accordingly an investigation was undertaken to evaluate this and 
other factors contributing to the phenomenon of chilliness, or clammi- 
ness, in wet fabrics. 
II. MATERIALS 


. 


Thirty-six samples of knit underwear fabric, prepared especially for 
experimental purposes, were used throughout ‘the investigation. The 
fabrics consisted of plain and special knits of various weights and 
thicknesses and of different fiber compositions, as described in tables 
1 and 2. Both unlaundered and twice-laundered samples were 
examined. 

TABLE 1.—Description of fabrics 














Nominal 
} | weight 
| of under- 
Sample Type of spinning | Type of knitting Fiber content of yarn clone 
| , 
| size 40) 
unlaun- 
dered 
| lb 
PE ee | Cotton_....--- ond Be craks sicaeupineete Be Emr flee 10% 
RR do : a CR do ls aliens a 75% cotton, 25% wool. wucaaws 10% 
Be Seates | French es Hane? “esereperrnr FO eg a See 10% 
 RREINRE SEE: FEES pee REE N” al e ar ee Zp , ..:.. san oun 10% 
RY ee” ae ch, BE. do. pa cea tans 25% cotton, 75% ES? } 10% 
F | nisendl | OEE ee ES Me is WOE. ck eens cate } 1034 
(Aid _.| Cotton... Sagas BOE a aes Peak meee en emapenn | 8 to 8% 
SSS SY eee eee a 75% cotton, 25% wool. .......----- 8 to 8, 
| AV £_ , SRSeeeeRee en | send SO ics letaebdo ican 50% cotton, 50% wool. -.....------ 8 to 8% 
ae Bee PEE eee wert do.........-....--.| 28% cotton, 75% wool. ..........-.- 8 to 84 
A ..do_. Mle: Atl, | Sinead Py NU cs As ow adenetnns ogee 8 to &&% 
| | Cotton...... vader acaeel ented | ESE Cees tg ESA Ge aed © apr eae near 5% to 6 
| ae Ls ASR ose Biivedd aitea athe Re: Dh old eetiadhlt 75% cotton, 25% wool....--.------ 5% to 6 
_ bo ee eae, Pas | Sa SLR alsa 50% cotton, 50% wool......--..---- 5% to 6 
2 ee eS eee pokinwel | teas ete Jantee’ 25%, cotton, 75% wool. ..-....--.-.-. 5% to 6 
eee } do Bei Ge eet SO WOE oek Seis st _ oe eae 5% to 6 
| | 
CA SEE. | Mesh stitch._...-..-- 50% ¢ on 50% wool. ..-.--- ane 8 to 8% 
CB... = — “| ae aee 58 CDs iawn vd dekeudaeessotebeess ..-| 5% to 6 
CC. _do pete .| Waffle stitch fe : o el i ge Co aria 10% 
i Filament_____- a a 100% acetate filament__........_- Light * 
<a ey ere gh Ce Shes ...| 100% viscose-rayon filament. ----.- Do. 
OP .i....! _, ES | ee pa Do. 
yee ---] EEE es ee Ys : PR Lae 50% « ane, 50% wool. ....-.-..-.- 10% 
PCRS ee” Se ac: acacia ae aS SEE “RARE NS CRE TRE eer 8 to 84 
Parco aie. ee ed SEP oon FF St GPRS ESA 5% to 6 
. ee | French.........--.- w-]-=2-20-nansooeeeeene-- 50% cotton, 50% wool.._.....----- 5% to 6 
tg, easel | ulMetce munis on Ed WERPRIREER FECIL, | SERRUMNSD TOA E: chy erasmeyaes 10% 
: SSRN | Cotton and French_- na Duofold type RE POST io pmo in, 100% wool out__.-- 104 
Sea ae wal .do 22-2--| 100% cotton out, mk fonncrg ete 10% 
SD. . inccd MUM cotta ndtcr cages | Interlock.............- b0% « cotton, 50% wool......-....-- 8 to 8% 
XE - me Rec Bee ie: “ited ee i GR ee 104 
| 
ae in. 1 i ra hal | Alternate rib feed.__..| 75% cotton, 25% wool_...-.--.---- 104 
> PRE citit? “PUR Ae SEE’. GMa eae 50% cotton, 80% Wool... /.....5-0- | 8 to 8% 
| an ent TS eR | Flat fleece lined. -.___- ee BS et Ee Ea Te 1044 
Acdsee 2 SAN ie Siege lee My WOE 2) dh wchion do nennenes 10% 
XJ .....--| OU S00. Rw 4S 100% cotton.......... RIS HERE 10% 
| } 











® The nominal weights of these fabrics in pounds per dozen undershirts of size 40 was not given. Their 
actual weights in ounces per square yard are as follows: CD, 3.4; CE, 4.2; CF, 5.4 


2 Figures in brackets indicate the literature references at the end of the paper. 
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TABLE 2.—Properties of twice-laundered fabrics 








| ; 
| Thickness Air perme 























| | 
| } 
Weight sshility » . 
Sample of iz fabrics t pressure | — 
fabrics pressure "op of | ransmission 
| of 0.1 Ib/in.| | 
. H20 | 
| | ozyd? | in. | Fi eminx fe) | Btu/(° haa “i 
tA 11.4 | 0. 063 1.2 
B Teen 12.3 | 073 | 82 124 
C SE 12.0 | 078 85 1. 20 
D Se Se 10.8 | . 066 126 1. 25 
B... .| 11.7 | .070 137 1. 24 
| } 
a 9.8 | 076 | 197 1.18 
7... 9.8 | 061 | 101 | 1. 30 
#7. 9.7 .072 | 109 | 1. 20 
. te Ree | 10.9 077 94 | 1. 20 
ja SATS | 9.5 | .073 195 | 1. 26 
ig. 0.9 | 067 215 | 1.20 
| Yat Oe 6.0 | . 048 180 | 1. 33 
| di A Tgp RS 6.3 | 055 | 193 1. 26 
@ {Brae 6.2 | 054 239 | 1, 25 
Q.. --| 6.4 | . 057 | 268 | 1. 24 
R.. 7.0 | 080 | 310 | 1.15 
CA. 8.3 068 | 238 | 1.32 
of DES | 4.9 | 055 300 | 1.32 
CC... 13.0 | . 108 144 | 1. 28 
a! Uae Ras i | 3.6 . 023 e249 | 1. 48 
| Oe meee 4,2 | . 028 © 242 1. 48 
s) gal SRS | 5.5 | . 032 © 106 1. 42 
RE 9.1 | . 060 102 | 1. 26 
aE se ae | 6.9 | . 044 195 | 1. 33 
FM.. g 5.7 | . 045 189 | 1.34 
| 
tennis 5.4 | . 040 56 | 1. 33 
pA RETRO STEN 9.9 | .070 164 | 1, 22 
ea oie’ 10.7 . 095 178 | 1.16 
So aS 10.9 091 162 | 41,07 
XD 10. 4 072 | 99 | 1. 24 
| 
XE 10.4 097 | 116 | 1.12 
XF 10.7 .077 | 106 | 1. 24 
ae << Sree 9.8 074 121 1. 29 
a. 405 8.9 115 81 | 1.04 
9.0 ‘089 204 1, 18 
| Mido tabs 8.2 075 | 504 | 1.31 








* Measured with the instrument described by Schiefer and Boyland [7]. 

» Measured with the instrument described by Cleveland [8]. 

¢ The values for the air permeability of samples CD, CE, and CF are for a pressure difference across the 
fabric of 0.05 inch of water, not 0.5 inch as used for all other samples. The air permeability of these samples 
was too high to measure at 0.5-inch pressure difference. 

4 With sample XC, a duofold fabric, the test was first carried out with the wool layer inside (i. e., against 
the heating plate of the apparatus) and then with the cotton layer inside. The results were 1.07 and 1.14, 
respectively. 


III. EXPERIMENTAL PROCEDURE 
1. SUBJECTIVE TESTS 


These tests were carried out in a room conditioned at 21.1° C and 
65-percent relative humidity. Samples of the fabrics, 3% inches square, 
were soaked in water at the temperature of the conditioning room, 
and then dried in a towel until they contained the desired amount of 
water. The fabrics were tested when containing 100, 150, and 200 
percent of moisture, based on their conditioned weights. All the fab- 
rics were used in subjective tests, although each fabric was not tested 
with every other fabric because of the large number of tests that would 
be required. They were tested two at a time, on five male subjects, 
by simultaneously placing one fabric on the right and the other on the 
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left arm (radial ulnar region of the forearm). The subject was thay 
asked whether he detected a difference in chilling effect between th: 
two arms, and if so, which fabric gave the greater sensation of chilling. 
Each fabric was left on the arm for 15 seconds, and then removed. 
because it was found that the subject’s immediate reaction appeared 
to be most reliable. Although some persons were able to perceive 
slight differences better than others, the subjective reactions of various 
individuals were in very good agreement. 


2. MEASUREMENT OF DROP IN SURFACE TEMPERATURE 


The apparatus used to determine the drop in surface temperature 
was essentially a constant-temperature oil bath (fig. 1). In order to 
provide conditions under which a substantial drop in “skin” tempera- 
ture could be measured, the heated surface was covered with several 
layers of asbestos sheeting and a layer of waterproof rubberized cloth, 
the total thickness of this cover being approximately %.inch. The skin 
temperature was measured by four copper-constantan thermocoupies, 
symmetrically placed on the outerside of the rubberized cloth, and 
held in position by a thin layer of cellulose nitrate cement. The 
thermocouples were connected in series, and the ice point (0° C) was 
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Fiaure 1.—Diagram of apparatus for measuring drop in “skin” temperature. 


U pper and left portions, schematicside yr: lower ge cross section. A, Outer box; B, thermoregulator; 
C, mechanical stirrer; D, heater; E, oil; F, rubberiz fabric; G, asbestos; H, metal container; I, wooden 
box; J, insulation; K, radiator; L, cheesecloth; M, fan; N, water tank; and O, pump, 
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ysed as the reference temperature. A switch arrangement made it 
possible to measure with a potentiometer each individual electromo- 
tive force, as well as the sum of the four individual values. The 
average Skin temperature was calculated from one-fourth of the value 
of the potential difference of the four units in series. To insure a 
highly constant skin temperature, the heating apparatus was placed 
in a small conditioning cabinet, which in turn was placed in a room 
maintained at 21.1° C and 65-percent relative humidity. Water at 
room temperature is circulated through the radiator, K; the fan, M, 
blows air at room temperature first over the radiator and then over 
the heating apparatus. Since the water in the radiator does not 
respond to the relatively short-period temperature fluctuations of the 
air in the room, the temperature within the cabinet is considerably 
more constant than that in the room. Using this device, the average 
temperature of the ‘“‘skin” surface was constant to +0.1° C. The 
internal temperature of the bath was held at 50.0°+0.05° C, and the 
average corresponding skin temperature was 37.5° C.* 

The experimental procedure consisted in first measuring the skin 
temperature when the apparatus, without a test fabric on it, was in 
equilibrium with its surroundings. Subsequently, a moist fabric was 
placed on the apparatus. The moist fabrics were prepared by soak- 
ing approximately 3.5- by 6-inch samples in distilled water, the tem- 
perature of which was the same as that of the conditioning room, and 
then pressing the samples gently between towels until they contained 
150 percent of moisture, based on the weight of the conditioned cloth. 
The moist fabric was placed on the apparatus so as to cover the ther- 
mocouples, and was covered immediately with a cheesecloth belt 
(approximately 4 in. wide) a. a at both ends (0.3 lb per end) 
in order to hold the sample in close contact with the heated surface. 
The change in surface temperature was followed as a function of time, 
the first reading being made 30 seconds after the initial contact of 
the moist fabric with the heated surface. Observations were made 
for 10 minutes and then discontinued, because it was found that fur- 
ther changes in temperature occurred only very slowly. 


3. “CONTACT” TESTS 


Since it is reasonable to expect that the chilling effect of wet fab- 
rics may be related to the extent of contact which they make with 
the skin, a method was devised to check this hypothesis. A moisture- 
sensitive paper was prepared by impregnating a soft paper‘ with a 
water-soluble dye of high tinctorial power,® according to the following 
procedure. A mixture consisting of 5 g of dye, 30 g of crude milled 
rubber, and 1 liter of Stoddard’s solvent was prepared by first dis- 
solving the rubber in the solvent, and then grinding together the dye 
and rubber solution in a mortar. To insure uniformity of distribu- 
tion, the mixture was subsequently placed on a mechanical shaker 
for half an hour. The dye is insoluble in the organic solvent and 
does not appear red under these conditions; the rubber aids in dis- 

4 Since a setup that would discriminate between fabrics with fairly similar thermal properties was required, 
the drop in skin temperature was accentuated by using a relatively hot skin (37.5° C) and a thick layer of 


insulation. The good correlation of these measurements with subjective tests demonstrates the validity of 
this procedure for the purpose of evaluating the fabrics. 
‘The paper was a soft sheet made from wood pulp (No. 1 alpha), containing a small amount of rosin 
Size (sizing value by the dry indicator method: 28 sec), basis weight 23 to 24 lb, 17 by 22 in.—500 sheets. 

' The dye was basic fuchsine (du Pont fuchsine concentrated powder, L7596). 
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persing the dye and imparts to the finished paper certain desirable 
properties. Strips of the paper were dipped in the impregnating 
mixture and hung up to dry at room temperature. The last traces 
of solvent were removed from the papers by placing them in a vacuum 
oven at 65° C for 15 to 20 minutes. After trimming the papers to 
the desired size, they were ready for use. 

The test for contact was carried out by placing pieces of the wet 
fabrics, 3% inches square, on the papers for 10 minutes, under a pres- 
sure of 0.5 lb per square inch. A piece of the sensitized paper was 
first laid on a dry glass surface. The wet fabric was then placed on 
top of the paper, and covered with a glass square of similar size, and 
weights to bring the pressure to the stated amount. After 10 minutes 
the fabrics and weights were removed. The papers, which were 
practically colorless originally, became red where the moist fabrics 
touched them, thereby indicating the extent of contact of the fabrics,’ 
Photometric measurements of the reflectance of the papers were sub- 
sequently made, thus making it possible to assign to each fabric a 
numerical value indicative of its contact at any given moisture con- 
tent. The measurements were made by means of a Martens pho- 
tometer with a filter to isolate the 560-millimicron line. 


IV. RESULTS AND DISCUSSION 


1. QUALITATIVE RELATION BETWEEN SUBJECTIVE AND OBJECTIVE 
TESTS 


This investigation was undertaken not to develop a standard tech- 
nique for evaluating the chilling properties of fabrics, but rather to 
develop a method for studying a phenomenon which has been ob- 
served empirically for a long time. Owing to the many complex 
factors which contribute to this phenomenon, an attempt will not be 
made to present the data on a quantitative basis. It is recognized, 
furthermore, that only one aspect of the larger problem of the thermal 
properties of fabrics has been studied. Nevertheless, several signifi- 
cant relationships were noted, and these will be discussed in more or 
less general terms, making no effort to account for small differences 
between the many and diverse fabrics included in the investigation. 

The results of the subjective tests were in very good agreement with 
the ‘direct measurements of drop in temperature and degree of con- 
tact. The results of a typical set of subjective experiments are given 
in table 3, and may be compared with the corresponding objective 
measurements of figures 2 to 21. When, as, for example, with sample 
A, the subject experienced pronounced chilling, a substantial drop in 
temperature was measured, and the fabric was found to make good 
contact with the skin. On the other hand, a smaller drop in tempera- 
ture was measured, and poor contact was observed in fabrics such as 
sample F, which did not produce a sensation of chilling to the subject. 
This correlation between subjective reaction and contact was especi- 
ally interesting when the fabrics were placed on the moisture-sensitive 
papers immediately after their removal from the arm, Invariably, 
when the subjects experienced pronounced differences in comfort, 
similar outstanding differences in contact were observed. Likewise, 


6 This test does not actually measure the contact which the fabrics make with the human skin, but for 
the purposes of the present investigation, the test is believed to be an adequate indirect measure of this 


property. 
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Figure 2.—‘‘Contact” impressions obtained with fabrics containing 200 
percent af moisture. 
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Figure 3.—‘‘Contact”’ impressions obtained with fabrics containing 200 
percent of moisture. 
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Figure 4.—‘‘Contact’’ impressions obtained with fabrics containing 200 
percent of moisture. 
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when the individuals perceived little or no difference in chilling effect, 
the sensitized papers showed that the two fabrics made similar con- 
tact. A typical series of contact papers, obtained when the fabrics 
contained 200 percent of moisture, is shown in figures 2 to 5. It is 
clear that, although all the fabrics contained equal percentages of 
water, some made better contact with the papers than others. Thus, 
for example, samples A, G, M, FM, amd XJ made very good contact, 
whereas Others, such as F, L, R, XC, and XJ, made poor contact. 
By increasing and decreasing the amounts of water in the fabrics 
their ‘‘tolerance’’ to moisture was estimated. It was found that even 
at 450 percent of moisture some of the fabrics (L, CB, XE, and XH) 
still did not completely color the papers. On the other hand, at 50 
percent of moisture vale fabrics A, B, G, M, CF, FM, and XJ colored 
the papers at all, and then only lightly. This variation in contact is 
shown graphically in figures 6 to 13, where the reflectance of the papers 
as measured photometrically is plotted against percentage of moisture 
in the fabrics. 

To facilitate a more detailed discussion, the fabrics may be grouped 
according to weight, construction, and fiber composition. 


TaBLE 3.—Subjective results obtained with fabrics containing 200 percent of moisture 
on the basis of their conditioned weights 








Samples com- 
pared 


Samples com- 


subiecti A 
pared Subjective reaction 


Subjective reaction | 
| 





Fand XI F somewhat colder. 





Qand FO 
Gand L.__....- 


i oe Som 
Aaet 8.6 i08508 


XFand XG... 





A much colder. 


.| M much colder. 


Neither produces much chill, but 
N slightly colder. 

Neither causes much chilling, 
but FO slightly colder. 

G much colder. 


Both rather cool, H slightly cool- 
er than K. 

FI slightly colder. 

A much colder. 

B colder. 

No difference noted. 

Both produce some chilling, but 
XF slightly colder. 








CH and CF ....- 
CE and CD___- 
XBand XC__.- 
FO and FI____.- 


FM and FC__-- 
XJ and XA_-.-_. 
CBand CA._--. 


CC and CD....-. 
XE and XH_-_-- 








CF slightly colder. 

Both cold, CD colder. 
XB cold, XC warm feel. 
No difference. 


FM colder. 

Both cold, but XJ colder. 

Very little difference, but CA 
may be slightly colder. 

CD much colder. 

Both comfortable, but XE very 
slightly colder. 
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FicurE 6.—Reflectance of ‘‘contact’’ impressions obtained with fabrics containing 


various percentages of moisture. 
10%-lb fabrics composed of various percentages of cotton and wool. 
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Figure 7.—Reflectance of ‘‘contact’’ impressions obtained with fabrics containing 
various percentages of moisture. 
& to 84-lb fabrics composed of various percentages of cotton and wool. 
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Fiaure 8.—Reflectance of “‘contact” impressions obtained with fabrics containiy 
various percentages of moisture. Fic 
5}4- to 6-lb fabrics composed of various percentages of cotton and wool. 





Thermal Properties of Moist Fabrics 239 





























70r 
CD-O-ACETATE 
CE-@-VISCOSE 
CF-@-NYLON 
60r 
4 
al 
< 
oO 
Ww 
eS 
uj 407 
« 
¢ 
— 
8 GF 
ys SOP 
i¢ ) 
20r 
0 
ior 
fe) ] | { ] 
100 200 300 400 


PERCENT MOISTURE 


Ficure 9.—Reflectance of ‘‘contact’”” impressions obtained with fabrics containing 
various percentages of moisture. 
Lightweight fabrics made of continuous filament acetate, viscose, and nylon. 
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Fiaure 10.—Reflectance of “‘contact’’ impressions obtained with fabrics containing Pr 


various percentages of moisture. 
1044-lb fabrics of different constructions, composed of 50 percent of cotton and 50 percent of wool. 8 
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9 Figure 11.—Reflectance of “contact” impressions obtained with fabrics containing 


various percentages of moisture. 
& to 84-lb fabrics of different constructions, composed of 50 percent of cotton and 50 percent of wool, 
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Figure 12.—Reflectance of ‘‘contact’’ impressions obtained with fabrics containing 
various percentages of moisture. 
1044-lb fabrics of different constructions and various fiber compositions. 





Thermal Properties of Moist Fabrics 








FM-@-FLAT— 100% COTTON 
FO-®-FLAT — 50% COTTON, 50% WOOL 


CB-©-MESH STITCH — 
50% COTTON, 50% WOOL 











) 


S 


rs} 
a 
= 
uj 40 
2 
© 
x 











i L r j 
100 200 300 400 
PERCENT MOISTURE 
MieurE 13.—Reflectance of “‘contact’” impressions obtained with fabrics containing 
various percentages of moisture. 
5¥- to 6-Ib fabrics of different constructions and various fiber compositions. 
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2. INFLUENCE OF FIBER CONTENT ON CHILLING EFFECT 


Figures 14 to 17 show the drop in temperature caused by laundered 
fabrics of comparable weight and construction but of different fibe 
compositions. The curves show a rapid initial drop in surface ten. 
perature’ during the first 30 seconds of the measurements, usually 
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Figure 14.—Drop in temperature upon application of moist fabrics to heated 
“skin’’ surface. 
10%-lb fabrics composed of various percentages of cotton and wool. 


followed by a slight increase in temperature, after which an approxi- 
mately steady state is reached. 
Since the thermocouple readings are influenced by the underside of 
the moist fabric as well as by the skin, the readings do not correspond 
recisely to the skin temperatures, This fact probably accounts in 
I rge measure for the unusual shape of the curves in the early portion 


1 In the course of measurements on the transmission of heat through textile fabrics, Rees [5] observed 8 
very small drop in temperature when fabrics conditioned at 25° C and 65-percent relative humidity were 
—_ in contact with a hot plate. The magnitude of the drop was greater for cotton than for woolen 
fabrics, 
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of the measurements. Thus, the rapid initial drop in temperature 
caused by the contact of the moist sample is followed by a rise in tem- 
perature, corresponding to an interval in which the heat supplied 
exceeds the heat lost by evaporation. When the loss of heat by evap- 
gration becomes sufficiently Eee: a decrease in the thermocouple 
temperature again ensues. he heavier fabrics, because of their 
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Figure 15.—Drop in temperature upon application of moist fabrics to heated 
“skin” surface. 
8 to 8%-lb fabrics composed of various percentages of cotton and wool. 





higher heat capacity and greater (absolute) content of water, cause 
alarger initial drop in skin temperature (compare figs. 14 to 17). On 
the other hand, a longer time is required for the evaporative losses 
to exceed the heat supplied in the case of the heavier samples. 

In general, the results given in figure 14 show a regular, progressive 
decrease in the chilling effect of the 10%-lb * fabrics with increase in 
their wool content. Similar results for samples of nominal weights of 
8 to 8% Ib and 5% to 6 lb are shown in figures 15 and 16, respectively. 
Sample R (fig. 16), which showed the smallest decrease in temperature 
ee 
‘The nominal weight of these fabrics is 1044 Ib per dozen suits of size 49. For brevity, they will hence- 


forth be designated *10}4-lb fabrics,” and a similar designation will be used for fabrics of other nominal 
Welghts. Although nominally the same, the fabrics differ appreciably ia actual weight, as shown in table 1. 
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of all the rib-knit samples, was characterized by a very fuzzy surface 
after laundering and was considerably thicker than the other fabrics of 
comparable weight. Comparison of the curves showing a drop iy 
surface temperature (figs. 14, 15, and 16) with those showing the re. 
flectance * of the contact papers (figs. 6, 7, and 8) reveals a good quali- 
tative relation between these two sets of measurements. 

Figures 9 and 17 show the reflectance of the contact papers and the 
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Ficure 16.—Drop in temperature upon application of moist fabrics to heated 


“skin’’ surface. 
5}4- to 6-lb fabrics composed of various percentages of cotton and wool. 








drop in temperature, respectively, of the lightweight, continuous 
filament underwear samples. These fabrics cause a big drop in tem- 
perature, but the differences between the acetate, viscose, and nylon 
fabrics are not large, and are of questionable significance. The difl- 


P * The more color produced on the papers as a result of contact with the wet fabrics, the lower is their re 
ectance, 4 
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wity in manipulating these very thin fabrics in the various tests may 
count in part for the low correlation of the contact and temperature 
measurements. 


3. INFLUENCE OF CONSTRUCTION ON CHILLING EFFECT 


Figures 18 and 19 show curves of the drop in temperature obtained 
sith fabrics of varied knit but of the same fiber content (50% cotton, 
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FiaurE 17.—Drop in temperature upon application of moist fabrics to heated 
“skin” surface. 
Lightweight fabrics made of continuous filament acetate, viscose, and nylon. 





50% wool). The curves of rib-knit fabrics (samples C and J) are 
included for comparison. Because of the complexity introduced by 
fabrics of different constructions, the apparent relation between drop 
in temperature (figs. 18 and 19) and extent of contact (figs. 10 and 11) 
is not so close as with fabrics of simple rib knit. However, fabrics 


XE and XC (terry stitch and duofold" with wool inside, respectively) 


“The duofold fabrics consist of two layers, one all cotton and the other all wool. 
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appear to be superior to the conventional rib knit. Since the duof| 
fabrics are identical, except that one has cotton and the other wool o 
the inside, a comparison of the results for these two samples clegy| 
shows the advantage of having wool next to the skin rather than ¢o} 
ton. Although fabrics which are superior at one moisture percentag 
are usually superior at another also, this is not always the case, as j 
shown by crossing of the reflectance curves. 
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Ficure 18.—Drop in temperature upon application of moist fabrics to heated “‘skin” 
surface. 
10-lb fabrics of different constructions, composed of 50 percent of cotton and 50 percent of wool. 





Figures 20 and 21 show the drop in temperature for two groups of 
fabrics of various constructions and fiber compositions, Figures 12 
and 13 indicate the contact made by the same fabrics. The results 
of the two sets of measurements show a good qualitative relation, 
the fabrics causing the least drop in temperature being the ones which 
make the poorest contact. Sample XJ, which is an all-wool fabric 
spun on the woolen system, is among the best of all the fabrics tested. 
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he marked improvement that a fleece lining makes in an all-cotton 
sbric is illustrated by the curve for XH. 

Throughout the investigation the moisture content of the samples 
sas calculated on a percentage basis. If the fabrics had been tested 
‘ith a given moisture content per unit area, the relative merits of the 
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Figure 19.—Drop in temperature upon application of moist fabrics to heated “‘skin’ 
surface 
§ to 8%4-Ib fabrics of different constructions, composed of 50 percent of cotton and 50 percent of wool. 





various fabrics may have been somewhat changed. It is believed, 
however, that a similar basic relation between subjective tests, direct 
Measurements of drop in temperature, and extent of contact, would 
still have been found. 





250 Journal of Research of the National Bureau of Standards 





°C 
(1) 


~“ 


w 
g 
ra 
é 
= 
E 
z 
3 
3 





XF-O- ALTERNATE RIB FEED— 
50% COTTON, 50% WOOL 


XJ- @ - RASCHEL— 100% COTTON 

XH-@- FLAT FLEECE-LINED — 
100% COTTON 

X1-6- RIB— 100% WOOL— 
WOOLEN SYSTEM 














j J 


10 
TIME IN MINUTES 


FicuRE 20.—Drop in temperature upon application of moist fabrics to healeo 
“skin” surface. 
104-Ib fabrics of different constructions and various fiber compositions. 
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FicurE 21.—Drop in temperature upon application of moist fabrics to heated 
“skin” surface. 


5\- to 6Ib fabrics of different constructions and various fiber compositions. 
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V. CONCLUSION 


Although it is difficult to rate on an absolute basis the “chillino 
effect’’ produced by the diverse fabrics which were included in these 
tests, several important conclusions can nevertheless be drawn— 
thus, a progressive improvement of the fabrics as their wool conten; 
is increased and the superiority of certain types of construction are 
clearly indicated. The extent of contact which the fabrics mak. 
with the skin appears to be the significant factor. The tests show 
that those fabrics which make the poorest contact cause the leas; 
chilling. From this point of view, the desirability of wool fibers which 
are highly crimped and possess long range elasticity, properties which 
permit a type of fabric construction which minimizes the extent of 
contact which can be made with the skin, is apparent. In contrast 
cotton which exhibits considerable plasticity when wet is less desirable 
from the same point of view. However, the results show that special 
types of construction, especially those which produce a napped, or 
fuzzy, surface reduce appreciably the contact which even wet cotton 
fabrics make and thereby lessen the chilling effect. It appears, there- 
fore, that in selecting fabrics for warmth, especially when these 
fabrics are to be worn next to the skin in cold climates and under con- 
ditions where considerable perspiration is produced, consideration 
a be given to the type of contact which they make with the 
body." 


The data in table 2 were made available through the kindness of 
W. D. Appel, Chief of the Textile Section, National Bureau of 
Standards. 
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